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oxide  and  associated  thermal  effects  along  with  heat  transfer  between  the 
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programmed  for  computer  solution  for  determination  of  ignition  limits  and  igniticn 
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discussed. 
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X#  INTRODUCTION 

Recent  years  have  witnessed1  a large  development  of  fuel-rich 
propellant  technology  to  produce  formulations  for  primary  rockets  in  air- 
breathing  propulsion  systems.  Two  types  of  systems  have  been  especially 
prominent,  air  augmented  rockets  and  external-burning  missiles  or  projectiles. 
In  an  air-augmented  rocket,  the  incompletely  combusted  products  of  a fuel-rich 
primary  motor  formulation  are  mixed  and  burned  with  ram  air  inside  a secondary 
chamber  at  a pressure  equal  to  or  less  than  the  air  inlet  recovery  pressure. 

In  external  burning,  the  fuel-rich  primary  motor  exhaust  is.  injected  directly 
into  the  free-stream  around  the  missile  and  burns  at  atmospheric  pressure 
(sea  level  or  altitude) . Eoth  types  of  systems  primarily  utilize  composite 
solid  propellants  with  ammonium  per- hlorate  oxidizer:  thus  non-me tali zed 

formulations  typically  produce  compounds  consisting  of  C,  H,  0,  N,  Cl,  and 
possibly  F,  while  metalized  systems  also  produce  B,  Al,  Mg,  or  combinations 
of  these  metals. 

An  essential  prerequisite  of  success  in  air-breathing  applications 
is  rapid  and  efficient  combustion  of  the  primal;,  exhaust  in  air.  It  is  of 
interest  to  spell  out  here  the  conditions  for  this  combustion.  In  ducted 
propulsion  the  pressures  in  the  secondary  range  from  15  or  20  atm  down  to 
perhaps  half  an  atmosphere  depending  on  the  flying  altitude  and  Mach  Number, 
and  the  residence  times  typically  are  a few  milliseconds.  In  external  burning 
the  pressures  range  from  1 to  0.1  atm,  and  the  residence  times  are  a fraction 
of  a millisecond.  Thus  it  is  easily  seen  that  these  air-breathing  systems 
often  place  stringent  demands  on  preignition  and  combustion  times  of  the 
primary  exhaust  fuel  products  with  air. 

Boron  is  a partic-larly  attractive  ingredient  foe  propellants  for 
air-breathing  rocket  applications  due  to  its  high  heating  value.  Generally, 
it  is  desirable  that  the  amount  of  boron  in  the  fuel-rich  propellant  be 
maximized,  subject  to  processability  and  primary  motor  ejection  efficiency 
constraints.  The  practical  upper  limit  of  boron  loading  is  approximately  50 
to  60  percent  by  weight.  At  these  loading  levels,  very  little  boron  is  burned 
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in  the  primary  motor  due  to  lack  of  sufficient  oxidizer  — thus,  most  of  the 
boron  must  be  .burnec  in  particulate  form  with  air  in  the  secondary  combustor 
for  its  heating  Value  to  be  utilized.  In  addition,  at  these  high  loadings, 
primary  motor  flame  temperatures  tend  to  be  low.  These  factors,  combined 
with  the  low  residence  times  available  for  mixing,  ignition,  and  combustion 
processes  in  the  secondary  chamber  have  led  to  afterburning  efficiency 
problems,  particularly  under  operating  conditions  which  result  in  afterburner 
pressures  of  50  psia  or  less. 

Substantial  prior  work  has  been  reported  on  ignition  and  combustion 
of  pure  boron  (1,2, 3, 4, 5).  However,  at  the  outset  of  our  project  there 
remained  a large  uncertainty  as  to  the  nature  (physical  state,  chemical  com- 
position) and  especially  as  to  combustion  characteristics  o(  the  condensed 
material  discharged  from  the  primary  motor,  which  is  certainly  not  pure  boron. 
The  primary  reaction  must  be  expected  to  coat  the  existing  particles  with 
reaction  products  and/or  to  generate  new  particles.  Agglomeration  of  par- 
ticles is  a definite  possibility. 

The  objective  of  the  present  work  has  been  to  characterize  the 
nature  of  the  primary  exhaust  from  fuel-rich  motors  containing  boron,  and  to 
study  its  ignition  and  combustion  characteristics,  especially  as  compared  to 
the  known  information  for  pure  boron.  Section  II  of  the  present  report 
describes  (a)  the  rocket-motor  firing  program  in  which  the  condensed  exhaust 
(CE)  from  two  different  propellants  was  collected,  and  (b)  the  results  of 
thermodynamic  computations  for  these  formulations.  Sections  III  and  IV 
contain  the  results  of  the  physical  (microscopy,  density  measurement)  and 
chemical  analyses  of  the  CE.  Section  V deals  with  combustion  of  the  CE. 
Finally,  Sections  VI  and  VII  report  the  results  of  modeling  studies  describing 
the  ignition  of  single  boron  particles  in  hot  gas  streams  and  propagation 
of  flames  in  boron-oxygen-nitrogen  dust  clouds. 


ROCKET  PROPELLANTS  AND  MOTORS 


II  . 

Two  propellants  were  studied.  Both  were  composite  ammonium  per- 
chlorate-binder formulations  containing  large  amounts  (almost  50%)  of  boron 
and  a few  percent  of  magnesium  and  aluminum.  ARCADENE  (Atlantic  Research 
Corporation  trade  name)  256A  is  a carboxy-terminated  polybutadiene  formula- 
tion and  ARCADENE  280  a polyester  formulation.  The  boron  used  in  both  formu- 
lations was  the  Kawecki  Company  material  with  nominal  average  particle  dia- 
meters of  3pm.  However,  the  diameters  range  widely  from  about  1 to  20pm  or 
more.  Thus  both  formulations  are  made  up  of  the  following  elements:  C,  H,. 

0,  N,  Cl,  B,  Mg,  Al.  In  addition  to  these,  the  256A  formulation  contains 
a small  amount  of  iron,  because  it  is  catalyzed  by  n-butyl  ferrocene. 

Two  motors  were  fired  with  each  formulation  into  a large  cylindri- 
cal tank  (5.2  m long, 2.9  m in  diameter)  which  was  first  evacuated  and  then 
filled  with  argon  to  0.35  atm  pressure  prior  to  each  run.  The  rocket  chamber 
pressure  and  the  tank  pressure  were  measured  continuously  throughout  each 
run. 

Data  pertaining  to  the  four  motor  tests  are  given  in  Table  I and 
Figure  1.  Each  propellant  grain  was  weighed  before  the  test;  as  was  the 
residue  remaining  inside  the  motor  after  the  test.  Comparison  of  the  third 
and  fourth  columns  of  Table  I in  which  the  respective  weights  are  listed, 
shows  that  the  ejection  efficiency  ranged  from  93.87.  in  Run  4 to  99.4%  in 
Run  2.  During  all  tests,  the  tank  pressure  rose  continuously  from  the  initial 
to  the  final  value.  The  low  final  value  recorded  in  Run  2 can  be  ascribed 
to  the  fact  that  the  burning  time  was  about  twice  that  of  other  runs,  thus 
allowing  more  cooling  and  condensation  during  the  run.  As  shown  in  Figure  1, 
each  propellant  was  tested  at  a high  and  a low  pressure.  In  three  out  of 
four  runs  the  chamber  pressure  remained  at  a constant  value  over  a substantial 
fraction  of  the  test  duration;  in  Ruu  1 it  kept  rising  over  most  of  the  burn- 
ing time,  possibly  due  to  clogging  of  the  nozzle. 

After  each  test  the  tank  filled  with  argon  was  kept  closed  for  about 
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1 hour  so  the  CE  would  settle  in  a large  (5.4  m surface  area)  metal  trough 
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at  the  bottom  of  the  tank.  The  trough  was  then  emptied  into  a bottle  by 
remote  action  without  opening  of  the  tank.  Thus,  the  CE  samples  were  not 
exposed  to  air  while  hot.  Thereafter,  ordinary  precautions  were  observed 
to  avoid  undue  exposure  to  air  and  humidity  (e.g. , the  samples  were  stored 
in  corked  vials  in  desiccators),  but  no  attempt  was  made  at  absolute  exclusion 
of  atmospheric  oxidants. 

A set  of  two  thermodynamic  equilibrium  computations,  "A"  and  "B", 
was  made  corresponding  to  each  experimental  motor  run.  The  results  are  given 
in  Table  II.  Chamber  pressures  in  runs  2,  3,  and  4 are  the  plateau  pressures 
of  Figure  1.  The  Run  1 chamber  pressure  of  10.2  atm  is  an  average  value  over 
the  duration  of  the  run.  The  exhaust  pressure  was  always  taken  to  be  0.5 
atm,  corresponding  to  an  average  experimental  value.  Computations  "A"  took 
into  account  most  species  which  can  reasonably  be  ecpected  to  appear  in  the 
products.  Computations  "B"  omitted  B^C  and  BN,  two  species  which  are  very 
important  thermodynamically,  but  which  are  sometimes  suspected  not  to  appear 
in  reality  for  kinetic  reasons.  The  species  used  in  the  computations  are 
listed  in  Table  III.  The  underlined  species  in  Table  III  are  arbitrarily 
defined  as  "major":  they  are  the  only  ones  predicted  to  occur  in  amounts 

exceeding  0.5  weight  percent  of  the  total  formulation  in  any  of  the  computer 
runs.  Hov/ever,  they  are  not  always  "major";  e.g.,  in  computations  "B",  Mg(g) 
does  not  appear  in  any  appreciable  amounts.  The  last  two  columns  in  Table  II, 
giving  percentages  of  condensable  products , include  not  only  liquids  and 
solids,  but  also  those  gases  which  condense  when  cooled  to  300°K,  all  in 
amounts  computed  in  equilibrium  under  exhaust  conditions.  In  other  words, 
the  exhaust  equilibrium  was  assumed  to  remain  frozen  during  mixing  with  argon 
and  cooling. 

The  most  abundant  gas  species  is  I^,  which  accounts  for  more  than 
50%  of  the  total  gas-  in  all  computations.  Two  other  permanent  gases  are  CO 
and  N2,  but  the  latter  is  absent,  or  virtually  so,  in  cases  which  "allow" 
formation  of  BN. 
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The  only  significant  liquid  species  under  exhaust  conditions  is 
I^O^.  In  computations  "A",  the  amounts  of  B2O.J  (1)  range  from  5 to  12  weight 
percent  of  the  total  formulation  (6  to  13%  of  CE) . In  computations  "B"  these 
amounts  are  substantially  higher,  10  to  18%  of  the  total  (12  to  21%  of  CE)  . 


The  major  difference  between  computations  "A"  and  "B"  concerns  the 
formation  of  solids.  When  full  equilibrium  is  assumed  ("A"),  most  of  the 
boron  in  the  exhaust  appears  as  solid  B^C,  and  a significant  fraction  as  solid 
BN.  The  amount  of  free  boron,  if  any,  depends  on  the  ratio  of  boron  to  carbon 
in  the  formulations:  virtually  all  of  the  carbon  is  found  in  C0(g)  and  B^C(s), 

so  if  there  is  excess  boron  left  (after  BN,  B^C,  B20^  and  CO  form),  a sub- 
stantial fraction  of  it  will  appear  as  B(s).  ARCADENE  280  has  relatively 
little  carbon  (high  B/C  ratio,  low  C/0  ratio),  so  computations  "A"  on  Runs  1 
and  3 show  some  B(s) , and  no  C(s).  Computations1^"  on  ARCADENE  256A,  which 
has  substantially  more  carbon  (low  B/C  ratio),  show  a small  amount  of  C(s), 
but  no  B(s).  In  computations  "B",  the  formation  of  boron  compounds  is  arti- 
fically  suppressed  to  a large  extent.  Thus  much  of  the  boron  (range  of  87 
to  92%  over  the  four  runs)  appears  as  B(s) , and  most  of  the  remainder  (7  to 
12%)  as  B20^(l);  as  mentioned  previously,  the  amount  of  B2°3^  is  substant- 
ially higher  in  computations  "B"  than  in  "A".  Carbon  in  computations  "B" 
appears  almost  entirely  aj  C's)  and  CO(g) , but  the  abundances  of  these  two 
species  are  very  different  in  the  two  formulations:  in  ARCADENE  256A,  more 
than  90%  of  carbon  appears  as  C(s).,  and  less  than  107.  as  CO;  in  ARCADENE  280, 
only  about  407.  appears  as  C(s)  and  607.  as  C0(g).  In  all  computations  vir- 
tually all  of  the  aluminum  appears  as  solid  MgA^O^,  which  is  always  a major 
solid  species  (about  7%  of  total  exhaust).  Magnesium  distributes  itself 
among  MgO,  MgA^O^  and  MgCl2(g). 

The  major  condensable  gases  predicted  by  the  computer  are  MgC^, 

^2^*2’  ®2^3’  HC1,  and  Mg.  The  amount  of  these  gases  appears  to  be 

mainly  a function  of  the  exhaust  temperature.  Since  temperatures  are  always 
higher  in  computations  "A",  these  computations  predict  more  gases  than  compu- 
tations "B".  The  highest  exhaust  temperature  is  that  of  Run  1,  computation 


"A";  the  condensable  gases  according  to  that  computation,  comprising  19%  of 
the  total  exhaust,  consist  mostly  of  boron-containing  species.  At  lower 
temperatures  (e.g.  computations  "B")  the  equilibrium  of  boron-containing 
species  shifts  toward  formation  of  substantially  more  B203(l) ; the  condensable 
gases  drop  to  roughly  10%  of  the  total  formulation,,  and  they  consist  mostly 
of  MgCl2  and  HC1. 


The  theoretical  characterization  of  the  CE  relevant  to  this  report 
can  be  summarized  as  follows:  (a)  11  of  the  boron  will  be  found  in  the  CE. 

The  fraction  of  the  CE  which  is  boron  (any  form)  is  about  50%  for  either  propellant, 
regardless  of  the  type  of  computation  ("A"  or  "B") . Quantitative  data 
on  the  amounts  of  boron  in  the  CE  of  each  run,  both  as  predicted  by  the 
thermodynamics  and  as  determined  experimentally,,  are  detailed  in  Section  IV. 

(b)  Most  of  the  boron  (roughly  90%  of  it)  is  found  in  three  solid  species: 

B,  B^C,  and  BN.  The  question  of  the  realistic  distribution  of  boron  among 
these  three, i.e.,  the  relative  merits  of  computations  "A"  and  "B",  must  be 
resolved  by  experiment,  (c)  The  remainder  of  the  boron  is  in  the  CE  in  the 
form  of  boron-oxygen-hydrogen  species.  While  some  of  these  are  rather  volatile, 
all  will  be  condensed  at  the  room  temperature,  (d)  The  remainder  of  the  CE 
consists  largely  of  MgAl^  and  MgCl2  in  computations  "A",  and  MgAl^,  MgCl2, 
and  C(s)  in  computations  "B". 
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III.  PHYSICAL  ANALYSIS  OF  THE  CONDENSED  EXHAUST 

1.  Microscopic  Examination 

The  condensed  exhaust  (CE) , collected  as  described  in  Section 
II,  was  examined  under  an  optical  microscope  with  magnifications  up  to  500x. 

The  untreated  CE  consists  of  very  large,  but  only  loosely  agglomerated  clusters. 
The  tumbling  of  the  clusters  with  soft  (plastic)  balls  breaks  up  the  mass  of 
the  CE  into  smaller  particles.  Microphotographs  of  CE  samples  (a)  as  col- 
lected and  (b)  after  tumbling  are  shown  in  Figure  2.  The  treated  material 
(Figure  2b)  consists  of  particles  ranging  from  about  3 to  50  pm  in  diameter. 
Larger  particles  in  that  range  are  no  doubt  agglomerates  which  could  be  broken 
up  by  continued  tumbling.  The  small  material  ,(<  10  um)  may  be  either  firm 
agglomerates  or  single  particles. 

Several  distinctive  features  can  be  seen  under  the  microscope. 
The  bulk  of  the  CE  is  dark  and  apparently  amorphous,  but  quite  a few  crystals 
are  also  visible.  The  crystals  are  either  white  or  colorless.  The  dark 
material  could  be  carbon  and  carboneous  material,  elemental  boron,  or  possibly 
B^C.  The  light  material  could  be  any  and  all  of  the  following:  B^O^  (c°l°r~ 

less),  B2°3'  ^0  (triclinic  crystals,  white),  MgCl2  (colorless  hexagonal), 
MgCl^*  H^O  (white),  MgAl^  (colorless),  and  BN  (white  hexagonal). 

2.  Density  Measurement 

We  determined  the  density  of  the  CE  samples  by  suspending  them 
in  inert  liquids  having  several  different  specific  gravities « The  liquids 
were  mixtures  of  bromoform  (sp.  gr.  2.89)  and  carbon  tetrachloride  (sp.  gr. 
1.595).  The  samples  were  found  to  consist  of  particles  having  different 
densities,  mostly  in  the  range  of  1.8  to  2.1  gm/cc. 

A complete  set  of  such  density  measurements  has  been  made 
on  the  CE  sample  collected  in  Run  1.  Photographs  of  several  suspensions 
are  shown  in  Figures  3,  4 and  5.  These  photographs  allow  interesting  con- 
clusions to  be  made. 
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Figure  3 shows  suspensions  of  the  CE,  as  collected  in  Run  1, 
in  liquids  having  specific  gravities  ranging  from  1.8  to  2„2.  It  can  be  seen 
that  all  particles  in  the  sample  have  specific  gravities  lower  than  2.2,  and 
some  even  lower  than  1.8.  Most  of  the  sample  appears  to  be  around  sp.  gr.  2.0. 

A CE  sample  from  Run  1 was  washed  in  a hot  mixture  of  methyl 
and  ethyl  alcohol,  a standard  procedure  for  removal  of  ^rora  boron  powder. 

Suspensions  of  the  washed  sample  are  shown  in  Figure.  4.  It  can  be  clearly 
seen  that  the  specific  gravity  of  most  of  the  washed  material  is  higher  than 
2.3.  It  should  be  noted  here  that  MgCl9,  a major  thermodynamically  predicted 
constituent  of  the  CE,  is  also  soluble  in  alcohol. 

For  the  purpose  of  comparison,  two  suspensions  made  with  the 
"off-the-shelf"  boron  powder  (Kawecki  Co.  material  used  in  all  our  propellants) 
are  shown  in  Figure  5:  virtually  all  of  the  powder  has  the  specific  gravity 

of  about  2.3,  i.e.  , distinctly  higher  than  the  unwashed  CE,  but  perhaps 
slightly  lower  than  the  washed  CE.  The  handbook  value  for  the  specific  gravity 
of  pure  boron  is  2.32-2.35.  Thus  the  2.3  value  for  the  commerical  powder 
is  very  reasonable,  especially  since  the  material  is  known  to  contain  about 
2 % B2°3*  wbich  has  a lower  specific  gravity. 

Our  conclusion  is  that,  in  the  CE,  low-density  substances  coat, 
or  are  firmly  attached  to,  solid  B (and  B^C,  if  any).  The  washing  removes 
these  substances,  and  leaves  a residue  consisting  largely  of  B (sp.  gr.  2.32  - 
2.35)  and  possibly  graphite  (sp.  gr.  2.26),  B^C  (sp.  gr.  2.54),  or  3N  (sp.gr. 
2.25).  The  possible  low-density  species  are,:  (sp.  gr.,  1.85), 

(sp.  gr.  1.49),  and  MgCl9*  6H,0  (sp.  gr.  1.56);  the  MgCl^  without  any  crystal- 
line water  has  the  specific  g' avity  of  2.33.  Since  carbon  is  insoluble  in 
alcohol;  the  fact  tnat  the  washed  CE  contains  little  low-density  material 
suggests  that  the  untreated  CE  contained  little  or  no  amorphous  carbon  (sp.  gr. 
1.8  to  2.1).  This  has  to  be  interpreted  in  the  light  of  the  thermodynamic 
computations:  the  full  computation  on  Run  1 (Section  II) , including  B^C  and 

BN  in  the  products,  predicts  no  solid  carbon;  if  B^C  and  BN  are  "suppressed," 
the  thermodynamic  computation  predicts  the  CE  to  contain  about  77>  carbon. 

This  comparison  suggests  that  computations  "A"  may  be  somewhat  more  realistic. 
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IV.  CHEMICAL  ANALYSIS  OF  THE  CONDENSED  EXHAUST 

The  CE  samples  from  all  four  motor  runs  were  analyzed  chemically 
for  boron.  The  analysis  was  performed  by  Ledoux  and  Company,  Teaneck,  N.  J. 
on  a subcontract  from  Atlantic  Research.  The  CE  samples  were  first  separated 
into  the  soluble  and  the  insoluble  portions  by  boiling  water,  and  then  each 
of  the  two  portions  was  analyzed  for  total  boron. 

The  results  of  the  chemical  analysis  are  summarized  in  Table  IV, 
which  lists  (a)  fraction  of  the  CE  which  is  boron  (any  form),  (b)  division  of  the 
CE  into  soluble  and  insoluble  portions,  and  (c)  fraction  of  the  total  boron 
in  each  of  these  portions.  In  each  case  we  have  listed  the  thermodynamically 
predicted  values  along  with  the  experimentally  determined  ones.  The  species 
assumed  insoluble  for  the  purpose  of  comparison  of  the  computation  with  the 
experiment  are  listed  in  Table  IV.  The  reader  should  consult  Table  III  to 
see  which  species  can  occur  in  any  appreciable  amounts  (>  0,5%  of  total  form- 
ulation). Inspection  of  Tables  III  and  IV  shows  that  the  insoluble  boron 
can  be  expected  to  be  found  in  B,  B^C  and  BN.  The  soluble  boron  may  come 
from  a number  of  boron-oxygen-hydrogen  species. 

Before  discussing  the  individual  columns  of  Table  IV  we  point  out 
that,  generally,  the  experimental  results  show  somewhat  larger  differences 
from  run  to  run  than  the  theoretical  predictions;  moreover,  the  deviation 
of  the  theory  from  the  experiment  is  not  always  in  the  same  direction  from 
run  to  run  (e.g.,  always  positive  or  always  negative),  but  generally  random. 

The  experimental  procedure  consisted  of  three  independent  steps:  propellant 

combustion  (including  ejection  efficiency)  , sample  collection,  and  chemical 
analysis.  Since  the  combination  of  all  three  contributes  to  the  final  results 
listed  in  Table  TV,  we  cannot  ascertain  the  reasons  for  the  random  deviation 
in  all  four  runs.  However,  the  combustion  of  Run  4 does  appear  to  have  been 
faulty.  Table  I shows  that  the  ejection  efficiency  of  that  run  was  unusually 
poor.  Table  IV  shows  both  an  exceptionally  low  amount  of  boron  in  the  CE 
and  an  exceptionally  low  amount  of  insoluble  matter  in  the  CE;  both  results 
indicate  that  in  Run  4,  for  some  unknown  reason,  a significant  fraction  of 
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the  boron  had  remained  in  the  motor  (the  residue  was  not  analyzed).  In 
view  of  these  atypical  results,  ie  shall  omit  Run  4 from  the  discussion  of 
trends  in  Table  IV. 

In  all  three  runs  under  comparison  (Runs  1,  2,  and  3),  the  measured 
amount  of  boron  in  the  CE  agrees  somewhat  better  with  computation  "A"  than 
with  computation  "B".  The  average  deviations  from  the  experiment  are  12%, 
for  computations  "A"  and  19%  for  computations  "B".  However,  it  must  be  pointed 
out  that  the  experimental  data  are  not  intermediate  between  the  two  computa- 
tions, as  might  he  expected  if  there  were  partial  formation  of  B^C  and  BN: 
rather,  both  theoretical  results  are  low  in  Run  2 and  high  in  Runs  1 and  3. 

The  partition  of  the  CE  into  soluble  and  insoluble  portions  is  pre- 
dicted distinctly  better  by  the  computation  "A"  in  Run  2,  slightly  better  by 
"B"  in  Run  1,  and  about  equally  well  in  Run  3.  Average  deviations  from  the 
experiment:  67«  for  computations  "A",  8%  for  computations  "B". 

The  percentage  of  insoluble  boron  is  predicted  better  by  "A"  in 
Runs  1 and  3,  and  equally  well  by  "A"  and  "B"  in  Run  2.  Average  deviations: 

12%  for  "A",  20%  for  "B". 

The  measured  amounts  of  boron  in  the  soluble  portion  of  the  CE  are 
predicted  quite  well  by  both  computations  in  Run  3,  but  very  poorly  in  Runs 
1 and  2.  Computation  "B"  appears  somewhat  better  of  the  two.  Average  dev- 
iations: 41%  for  "A",  29%  for  "B". 

We  conclude  that  according  to  both  the  experimental  data  and  either 
one  of  the  two  computations  only  about  50%  of  the  total  CE  is  boron.  Beyond 
that,  one  can  see  that  three  out  of  four  questions  which  can  be  answered  by 
the  results  of  the  chemical  analysis  - namely,  the  amount  of  boron  in  the  CE, 
the  fraction  of  the  CE  which  is  insoluble,  and  the  amount  of  boron  in  the 
'.nsoluble  portion  - are  predicted  reasonably  well  (order  of  10%)  by  the  full 
thermodynamic  computation  (column  "A"  in  Table  IV) . Computations  "B"  are 
somewhat  less  successful,  but  it  appears  that  they  also  give  a rough  idea  of 
the  composition  of  the  CE.  Indeed,  it  is  possible  that  the  assumptions 
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inherent  in  the  theoretical  comparisons.  - e.  gv,  the  freezing  of  the  exhaust 
equilibrium  - introduce  more  of  an  uncertainty  than  the  difference  between 
the  two  sets  of  computations.  On?.y  an  extensive  chemical  analysis  of  several 
species  (Mg,  B^C,  C,  and  others)  could  resolve  the  problem  of  the  detailed 
composition  of  the  CE. 
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V.  COMBUSTION  OF  THE  CONDENSED  EXHAUST 

As  described  in  Section  III,  and  shown  in  Figure  2,  the  ball-milled 
CE  consists  of  many  small  particles  with  sizable  agglomerates  interspersed. 
Prolonged  ball-milling  and  sieving  through  a 20  pm  mesh  screen  produced  a 
fine  powder  consisting  of  particles  having  a fairly  uniform  size  distribution 
with  average  diameters  of  3 to  4 pm.  By  careful  screening  between  30  and 
35  yin  we  also  found  it  possible  to  separate  coarse  powder  samples  consisting 
of  firm  agglomerates  which  could  be  treated  as  single  particles. 

Combustion  of  the  coarse  CE  samples  was  studied  by  a single-particle 
technique,  utilizing  a gas  burner,  which  was  previously  developed  and  exten- 
sively used  at  Atlantic  Research  for  combustion  studies  of  aluminum  (6,  7,  8), 
beryllium  (9 , 10)  and  boron  (2).  Since  particles  smaller  than  10  or  20  pm 
_in  diameter  are  difficult  to  handle  and  observe  singly,  we  adapted  the  gas- 
burner  technique  for  a study  of  particle  clouds.  The  newly  developed  tech- 
nique was  used  to  study  ignition  and  sombustion  of  the  fine  CE  powder. 

A.  The  Gas  Burner 

The  gas-burner  used  for  both  the  work  with  single  particles, 

30  to  40  pm  in  diameter,  and  clouds  of  fine  particles,  3 to  4 pm  in  diameter, 
is  shown  in  Figure  6 . In  the  single-particle  work,  the  powder  is  introduced 
directly  into-  the  f location  chamber;  the  lower  chamber  is  not  used. 

This  apparatus  was  described  in  detail  in  Ref.  6.  The  particles  are  entrained 
through  the  hypodermic  needle  by  an  inert  carrier  gas,  usually  helium,  and 
thus  introduced  rapidly  into  the  stream  oi  a burner  gas  having  known  composition 
and  temperature.  The  linear  velocity  of  the  burner  gas  is  of  the  order  of 
1 cm/msec i The  effect  of  the  inert  gas  on  combustion  of  particles  can  be 
neglected,,  because  the  hypodermic  tube  is  only  0.25  mm  in  diameter.,  and 
-ignition  of  particles  takes  place  several  centimeters  above  the  burner  plate. 

The  modification  of  the  burner  to  allow  combustion  studies  of 
clouds  consisting  of  small  particles  has  been  accomplished  by  (a)  the  addition 
of  the  lower  chamber  -t  and  (b)  enlargement  of  the  hypodermic  tube  to  0.6  mm 


in  diameter.  The  purpose  of  the  lower  chamber  is  to  break  up  particle  clusters, 
or  to  prevent  them  from  forming  in  the  first  place.  This  is  done  by  the  intro- 
duction of  a high-velocity  tangential  jet  of  the  carrier  gas  which  stirs  up 
the  powder  and  the  plastic  balls  with  which  the  chamber  is  packed.  In  this 
way  a steady  stream  of  powder  can  be  introduced  into  burner  gases  for  pro- 
longed periods  of  cime.  Our  experiments  usually  lasted  5 to  15  minutes. 

One  disadvantage  of  the  continuous  cloud  combustion  experiments  is 
that  the  effect  of  the  carrier  gas  on  ignition  of  particles  cannot  be  neglected. 
This  is  so  primarily  because  the  preignition  delays  of  small  particles  are 
very  short.  As  will  be  described  later  in  this  report,  the  CE  particles  ignite 
when  their  temperature  approaches  2000°K.  Heat  transfer  estimates  show  that 
heatup  times  from  room  temperature  to  2000°K  of  3.5  jum  diameter  -particles-,, 
injected  into  gases  having  temperatures  of  2200  to  2500°K,are  only  of  the 
order  of  10  to  100  psec.  The  laminar  diffusion  mixing  time  for  a 6 mm  dia- 
meter jet  will  be  of  the  order  of  0.5  msec.  Thus,  the  preignifcion  delays  and 
the  carrier-gas  mixing  times  overlap,  so  a quantitative  study  of  ignition 
delays  would  require  detailed  heat  transfer  and  mixing  calculations.  Out 
work  did  not  include  such  calculations.  On  the  other  hand,  studies  of  .post- 
ignition burning  times  are  not  seriously  affected  by  the  carrier  gas  mixing 
process,  because  the  temperature  of  the  gas  adjacent  to  the  particle  must 
approach  the  burner-gas  temperature  (i.e.,  the  mixing  must  be  nearly  completed) 
before  ignition.  Moreover,  in  the  actual  experiments  described  in  Section 
V-C,  the  carrier  gas  was  air,  and  the  mole  fraction  of  oxygen  in  the  burner 
gases  was  X = 0.2  in  most  experiments.  Thus,  the  mole  fractions  of  C^,  the 
most  important  oxidant,  were  matched  in  the  two  gases.  Since  the  carrier 
air  was  initially  cold,  the  effective  gas  temperature  during  combustion  was 
not  known  accurately  (it  was  probably  a little  below  the  gas-burner  temper- 
ature), but  the  gas  temperature  is  not  a dominant  parameter  during  self-sus- 
tained combustion;  the  mole  fraction  of  the  oxida»t  is.  Thus,  the  experiment 
is  well  controlled  as  far  as  post-ignitic-n  processes  are  concerned. 
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B,  Single-Particle  Experiments 

The  single-particle  technique  has  been  used  for  a study  of 
both  the  ignition  and  the  subsequent  combustion  of  a 30-40  pm  CE  sample  col- 
lected in  Run  4.  The  gas  burner  was  -run  on  propane-oxygen- nitrogen  mixtures. 

The  burner  gases  consisted  of  1.4  - 207,  0^,  most  observations  having  been 
made  with  about  207.  O2,  15  - 207.  F^O,  and  10  - 127.  CO^;  the  remainder  was 
mostly  N^.  Gas  temperatures  were  varied  from  3850  to  2300°K.  Velocities 
of  burner  gases  were  about  1G00  cm/sec,  and  those  of  the  carrier  (helium) 
about  100  cm/sec. 

Quantitative  ignitior  and  combustion  parameters  of  neat  boron 
powder  in  the  3C'  to  40  pm  range  are  well  known  from  our  previously  published 
work  (2).  Since  it  is  of  obvious  interest  to  compare  the  combustion  character- 
istics of  the  CE  to  those  of  boron,  we  shall  now  briefly  state  the  relevant 
results  of  Ref.  2.  The  entire  history  of  the  boron  particles  (d  = 35  pm) 
in  the  burner  gases  (2240°K,  197.  O2,  167.  H„0)  can  be  divided  into  three  phases: 
a dark  pre-ignition  period  (ca.  8 msec) , followed  by  a peculiar  two  stage 
combustion,  the  durations  of  the  two  periods  being  t^  = 4 rnsr:  and  t2  = 11  msec. 
The  first  combustion  stage  begins  when  the  particle  temperature  is  about 
1850°K.  The  first  stage  is  generally  believed  to  correspond  to  slow  combustion, 
impeded  by  a layer  of  2iquid  ^O^.  Most  of  the  mass  of  boron  is  consumed 
during  the  rapid  second  stage*  controlled  by  the  gas-phase  diffusion  of  oxi- 
dants toward  the  particle;  thus  t2  approximates  the  total  diffusion-combustion 
time  of  the  particle,  t2  ~ t^.  The  minimum  temperature  of  the  burner  gases 
necessary  for  ignition  of  particles  was  found  to  be  1880  ± 20°K. 

The  CE  agglomerates,  30  pm  < d < 40  pm,  were  found  to  burn 

as  single  particles:  there  is  no  separation  of  particle  trajectories.  Thus, 

we  can  make  a quantitative  comparison  of  the  CE  combustion  to  combustion  of 

boron.  When  the  CE,  30  pm  < d < 40  pm,  is  burned  in  the  same  burner  gases 

as  the  35  pm  boron  (T  = 2240°K,  197,  0„,  167.  H_0  - see  previous  paragraph), 

S 1 i. 

the  following  differences  are  observed: 
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1.  The  CE  "particles”  give  well-defined  burning  times,  but 
there  is  no  sign  of  two-stage  combustion. 

2.  The  average  burning  time  of  the  CE  is  t^  = 10  i 1.5  msec. 

3.  The  CE  ignites  substantially  lower  in  the  burner  gases  than 

boron:  we  estimate  that  the  preignition  delay  of  the  CE  is 

30  to  507.  lower  than  for  boron  particles. 

In  addition  to  the  tests  at  constant  gas  temperature  and  con- 
stant mole  fraction  of  0 we  ran  two  series  of  experiments  in  which  these 
two  parameters  were  varied.  Variation  of  temperature  at  constant  mole  fraction 
of  0^  (ca.  0.2)  showed  that  the  minimum  gas  temperature  necessary  for  ignition 
of  the  GE  is  1850  ± 25°K.  Variation  of  0o  content  at  the  approximately  con- 
stant  gas  temperature  of  2200  to  2300°K  shows  that  the  CE  continues  to  burn, 
although  very  slowly,  even  at  the  lowest  mole  fraction  of  0^  which  we  tried, 
0.014. 

Since  the  burning  time  of  a particle  is  a strong  function  of 
its  diameter,  and  since  we  have  no  accurate  value  for  the  average  diameter 
of  the  CE,  the  two  values,  t^  = 10  msec  for  the  CE  and  t^  ~ = 11  msec  for 

boron,  must  be  considered  the  same  within  the  experimental  error. 

C.  Cloud  Combustion  Experiments 

The  cloud  technique  for  combustion  of  small  particles,  des- 
cribed in  Section  V-A,  yields  average  preignition  delays  t.  and  burning  times, 
t^.  A cloud  of  burning  CE  particles  (Run  1)  is  shown  in  Figure  7.  Inspec- 
tion of  this  photograph  will  show  that  both  t^  and  are  approximate,  because 
in  a burning  cloud  one  cannot  define  clear  points  of  onset  and  cessation  of 
combustion.  Furthermore,  since  the  preignitior.  delay  necessarily  overlaps 
the  carrier-gas  mixing  time  (see  Section  V-A  for  discussion) , one  cannot 
obtain  cloud  ignition  temperatures  from  measured  preignition  delays  as  in 
single-particle  work;  one  can  only  determine  relative  ignitabilities  of  two 
clouds  by  comparing  the  preignition  distances  in  two  replicate  experiments. 

Two  CE  samples,  from  Runs  1 and  4,  have  been  studied  by  this 
technique.  The  results  obtained  with  these  two  samples  are  indistinguishable 
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from  each  other.  In  the  subsequent  sections  the  "CE  sample"  means  either 
Run  1,  or  Run  4,  or  both. 

1.  Particle  Sizes 

The  CE  samples  used  in  the  cloud  experiments  consisted 
of  the  material  which  had  passed  through  a 20  pm  mesh  screen.  However,  since 
the  gas  burner  operates  in  such  a way  that  the  material  is  subjected  to  pro- 
longed swirling  with  plastic  balls  (see  Section  V-A  and  Fig.  6),  there  is 
no  assurance  that  the  particle  size  distribution  injected  into  the  burner 
gases  is  the  same  as  in  the  original  (sieved)  sample.  We  therefore  sampled 
the  unburned  material  after  injection  into  the  burner  gases,  i.e.,  at  a point 
between  the  injection  orifice  but  below  the  point  of  ignition.  Microscopic 
examination  of  the  sampled  material  revealed  that  the  particle  diameters 
were  in  the  range  between  2 and  7 pm.  The  average  diameter  was  about  3.5  pm. 

In  addition  to  the  CE,  two  other  samples  were  studied  by 
the  cloud -combustion  technique  for  the  purpose  of  comparison:  an  amorphous 

boron  sample,  supplied  by  U.  S.  Borax,  Inc.,  stated  to  be  987.  pure,  and  the 
material  supplied  by  Kawecki,  Inc.  The  Borax  powder  particle  diameters  are 
in  the  narrow  range  of  1.0  i 0.2  pm.  The  Kawecki  powder  has  a wide  distri- 
bution of  particle  diameters,  ca.  1 to  20  pm,  with  the  number  average  of 
3 pm.  The  Kawecki  boron  is  used  both  in  ARCADENE  256A  and  280  propellants. 

2.  Ignition 

A technique  for  direct  measurement  of  minimum  gas  temper- 
atures necessary  for  ignition  of  particles  had  been  developed  previously  in 
connection  with  our  single-particle  work  (2):  a stream  of  particles  is 

continuously  injected  into  burner  gases,  the  temperature  of  which  is  gradually 
decreased  until  the  particles  cease  to  ignite.  This  technique  has  now  been 
successfully  applied  in  the  cloud-combustion  work.  The  mole  fraction  of  0^ 
in  these  ignition-limit  experiments  was  kept  about  0.2. 

The  limiting  temperature  was  found  to  be  the  same,  within 
the  experimental  error,  for  the  Borax  sample  and  for  the  CE:  1980  i 20 °K. 

It  will  be  noted  that  this  is  100  - 150°  higher  than  the  analogous  temperature 
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limit  for  30  - 40  pm  particles  of  either  the  CE  or  crystalline  boron  (Section 
V-B  of  this  report  and  Ref.  2) . 

Measurements  of  the  minimum  gas  temperatures  necessary  for 
ignition  of  powders,  which  we  shall  refer  to  as  "temperature  limits"  from 
now  on,  have  resulted  in  two  conclusions:  (a)  temperature  limits  for  clouds 

of  small  (<  10  urn)  particles  are  100  to  150°K  higher  than  for  large  (30  to 
40  urn)  particles;  (b)  regardless  of  the  particle  size,  temperature  limits 
are  the  same  for  the  CE  and  for  boron. 


We  shall  first  discuss  the  second  conclusion:  similarity 

of  the  temperature  limits.  This  result  is  surprising  for  two  reasons.  First, 
as  discussed  in  Section  IV,  only  about  half  of  the  CE  is  boron.  The  remainder 
consists  of  substances  which  have  different  ignition  temperatures.  Referring 
to  Tables  II  through  IV,  we  can  see  that  ignition  temperatures  of  some  expected 
ingredients  will  be  higher  than  those  of  boron  (B^C,  BN,  Al) , while  others 
will  be  lower  (Mg,  MgC^,  Fe,  C)  . One  would  expect  a priori  that  ignition 
of  any  appreciable  fraction  of  the  CE  would  lead  to  ignition  of  the  entire 
sample.  The  experimentally  determined  ignition  temperature  of  the  CE  perhaps 
suggests  that  the  amounts  of  ignitable  ingredients  are  relatively  small.  We 
have  no  chemical  analysis  of  the  CE  ingredients  other  than  boron,  nor  do  we 
know  the  critical  degree  of  admixture  of  ignitable  materials  required  to  lower 
the  ignition  temperature  of  samples  below  that  of  boron. 


The  measured  ignition  temperature  of  the  CE  also  appears  sur- 
prising at  the  first  glance,  because  the  observed  preignition  delays  are 
substantially  lower  than  for  boron.  This  result,  however,  can  be  rationalized. 
First,  we  point  out  that  the  density  of  the  CE  is  lower  than  that  of  boron 
(Section  III) . Now  the  heating  rate  of  a particle  is  inversely  proportional 
to  its  density: 
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where  p is  density,  C heat  capacity,  and  1c  thermal  conductivity;  the  subscripts 
p and  g refer  to  particle  and  gas  respectively.  This  means  that  the  heat-up 
time  of  a particle  to  any  specified  temperature  will  be  directly  proportional 
to  its  density.  Since  the  measured  densities  of  the  C£  are  10  to  20%  lower 
than  that  of  boron,  the  heat-up  times  will  be  lower  by  the  same  factor. 

In  addition  to  the  density  effect,  there  may  be  a more  subtle 
difference  in  the  pre-ignition  processes  between  the  CE  and  boron.  The  CE 
particles  are  almost  certainly  agglomerates,  so  they  are  apt  to  have  much 
larger  surface  areas  than  crystalline  boron,  which  was  studied  in  Ref.  2. 

Since  during  the  pre-ignition  stage  the  rate  process  will  usually  not  be 
controlled  by  gas -phase  diffusion  toward  the  particle,  but  rather  by  the 
global  surface  kinetics  - including  the  processes  of  diffusion  across  a sur- 
face layer,  vaporization,  and  chemical  reaction  (11.  12)  - the  pre-ignition 
rate  will  be  proportional  to  the  effective  surface  area.  The  large  surface 
area  of  the  CE  may  be  expe_ted  to  have  two  observable  effects:  (a)  the 

ignition  "runaway"  temperature,  and  therefore  the  pre-ignition  delay  will  be 
decreased;  and  (b)  the  duration  of  the  first  combustion  stage.  observed 
in  the  case  of  boron)  will  also  be  decreased.  As  described  earlier  in  this 
section,  a substantial  decrease  in  the  pre-ignition  delay  of  the  CE  as  com- 
pared to  boron  was  indeed  observed.  Concerning  the  second  expected  effect, 
decrease  in  t^,  we  have  observed  no  first  stage  at  all  in  the  case  of  the  CE. 
Since  a decrease  in  t^  by  a factor  of  about  4 would  suffice  to  make  the  first 
stage  elude  the  observation,  a large  increase  of  the  effective  surface  may 
account  for  the  effect;  however,  this  is  a matter  of  conjecture. 

We  now  turn  briefly  to  the  other  quantitative  conclusion 
regarding  the  ignition,  namely  that  temperature  limits  of  clouds  consisting 
of  small  particles  are  higher  than  for  large  single  particles.  A detailed 
theoretical  analysis  of  the  boron  ignition  (11,  12,  14)  shows  that  the  process 
is  complex  involving  problems  of  heat  transfer,  gas -phase  and  surface-layer 
diffusion,  and  vaporization.  Examination  of  Figure  17  in  Section  VI  (a 
reprint  of  Reference  12)  indicates  that  2 micron  particles  are  predicted  to 
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have  a minimum  required  surroundings  temperature  for  ignition  approximately 
90°K  higher  than  that  predicted  for  30  micron  particles,  consistent  with  the 
experimental  observations.  (This  is  not  conclusive,  however,  since  the  pre- 
dictions are  in  both  cases  for  single  particles  and  t’^e  cloud  effect  on 
minimum  required  surroundings  temperature  for  ignition  is  not  well-defined.) 

3.  Burning  Times 

The  cloud  combustion  technique  has  been  used  for  burning 
time  measurements  of  the  CE,  cf  the  Borax  boron,  and  of  the  Kawecki  boron. 

The  burning  time  t^  was  taken  to  be  equal  to  the  length  of  the  burning  cloud 
(see  Fig.  7),  divided  by  the  linear  velocity  of  the  burner  gases.  Since  the 
luminosity  of  the  cloud  both  develops  and  fades  out  gradually,  an  approximate 
visual  judgement  concerning  the  length  of  combustion  must  be  made.  The  vel- 
ocity of  the  burner  gases  can  be  calculated  from  known  (metered)  input  gas 
flow  and  known  flame  temperature;  it  can  also  be  obtained  by  measurement  of 
single-particle  velocities.  The  calculated  and  the  measured  values  usually 
agree  within,  about  25%.  We  estimate,  that  the  combined  uncertainties  of  the 
cloud  length  and  particle  velocity  may  conceivably  lead  to  an  error  in  deter- 
mination or  t^  of  up  to  a factor  of  two,  but  most  values  are  probably  better 
than  that. 

The  average  results  of  three  cloud-combustion  measurements 
are  shown  in  Figure  8.  Hie  burner  gas  properties  in  all  three  experiments 
were:  T = 2240 °K,  X(O^)  = 0.19,  X^O)  =0.16.  A single  point  for  the  t^ 

of  a spherical  boron  particle,  d = 37  um,  is  given  in  the  same  plot.  Tile 
single-particle  point  is  a small  extrapolation  of  an  accurate  measurement 
(see  Tabic  VI  of  Ref,  5)  in  a somewhat  different  burner  gas:  T = 2280°K, 

X(02)  = 0.23,  X(H20)  =0.  Even  with  the  extrapolation,  the  values  of  both 
t^  and  of  d are  no  doubt  much  more  accurate  for  the  single-particle  point 
than  for  the  cloud  measurements.  In  view  of  the  uncertainties  in  the  cloud 
measurements,  one  should  be  careful  about  assigning  a specific  numerical  value 
to  the  slope  n in  Figure  8.  However,  an  inspection  of  the  figure  shows  that 
even  if  the  cloud  data  are  in  error  by  a factor  of  2,  n will  still  be  closer 
to  unity  than  to  a value  of  2,  predicted  by  single-particle  diffusion  theories 
(13,  14). 
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In  addition  to  experiments  run  with  burner  gases  having 
a constant  relatively  high  mole  fraction  of  oxygen,  XCC^)  = 0.2,  cloud  experi- 
ments were  made  with  decreasing  values  of  X(O^)  both  with  the  CE  and  with 
the  Borax  boron  powder.  Gas  temperatures  in  all  the  tests  were  kept  at  about 
2300°K.  Only  visual  observations  were  made.  It  was  found  that  both  powders 
burned  about  equally  vigorously  down  to  X CO^)  = 0.14.  The  combustion  bright- 
ness decreased  rather  sharply  between  X(0,^)  =0.14  and  0.12,  and  then  very 
gradually  down  to  about  ^(O^)  = 0.05  with  both  powders.  At  very  low  mole 
fractions  of  oxygen,  0.015  to  0.05,  the  CE  powder  burned  very  weakly,  and 
the  boron  powder  hardly  at  all. 

The  single-particle  experiments  indicate  that  the  burning  times 
t^  of  the  CE  are  about  the  same  as  for  boron.  This  appears  a priori  reasonable, 
because  one  would  expect  the  limiting  process  to  be  the  gas-phase  diffusion 
of  oxidants  regardless  of  the  chemical  composition.  .Higher  volatility  of 
some  expected  CE  ingredients  should  accelerate  the  burning  rate  of  particles. 
Thus  if  the  CE  is  rich  in  species  such  as  Mg,  MgCl^,  and  Fe,  one  should  expect 
it  to  burn  somewhat  faster  than  boron.  On  the  other  hand,  if  it  contains 
large  amounts  of  B^C  (and  possibly  BN),  it  may  burn  more  slowly:  B^C  is  known 

to  burn  more  slowly  than  boron  (15)  , and  BN  has  been  observed  to  ignite  only 
with  difficulty.  In  either  case,  the  change  in  t^  of  boron  caused  by  adulter- 
ation is  probably  not  sufficient  to  be  clearly  demonstrated  by  our  experiments. 
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Boron  Particle  Ignition  in  Hot  Gas  Streams 

MERILL  K.  KING 

Atlantic  Research  Corporation , Alexandria.  Virginia 

Abstract — A model  of  boron  ignition  treating  the  inhibiting  effect  of  a boric  oxide  coating  has  been  developed. 
Transient  differential  equations  describing  the  generation  and  removal  of  the  oxide  and  associated  thermal  effects 
along  with  heat  transfer  between  the  particle  and  surroundings  have  been  derived,  converted  to  difference  form 
and  programmed  for  computer  solution  for  determination  of  ignition  limits  and  ignition  times.  Predictions  for 
particles  studied  by  Macek  in  a flat-flame  burner  have  been  compared  with  experimental  results-  agreement  is  good 


. irroundings  radiation  temperature  < 
by  a stability  analysis  to  determine  ignition  limits — excellent  agreement  is  found  between  results  obtained  with 
the  stability  analysis  and  the  numerical  analysis. 


NOMENCLATURE 

A,  B\  C,  D, 

constants  in  stability  analysis 

Qitx 

E,F 

dependent  only  on  initial  particle 

sire  and  such  boundary  conditions  as 
pressure,  oxygen  mole  fraction,  ana 

QitXi 

ambient  temperature 

R 

cl>lll 

liquid  boron  heat  capacity  (cal/gm  °K 

R' 

cl'U.S 

solid  boron  heat  capacity  (cal/gm  I) 

cfltjOj 

iiquid  boron  oxide  heat  capacity 

fin 

(cal/gm  °K) 

/ 

fraction  of  boron  in  the  liquid  phase 

It 

gas-pat  tide  heat  transfer  coefficient 
(cal, ’em2  sec  °K) 

Hr: 

k 

mass  transfer  coefficient  for  transport 
of  boric  oxide  gas  from  particle  to 

fin 

frec-slrcam  (gm-mol/cm2  atm  sec) 

0* 

M,  A’,  S,  U 

constants  defined  byequations  (20-23) 
fixed  by  A.  B,  C,  D,  E.  F.  T„  and  A' 

RGEN 

(.1/IF).. 

boron  atomic  weight  (gm/gm-mo!) 

REVAP 

(■UllOu.o. 

boric  oxide  molecular  weight  (gm' 

gm-mol) 

T, 

Nu 

Nussell  Number 

Tr,. 

P 

total  pressure  (atm) 

Tr 

oxygen  partial  pressure  in  free  stream 
(atm) 

^lljOj.-urf 

boric  oxide  partial  pressure  adjacent 

Ts 

to  particle  surface  (atm) 

boric  oxide  vapor  pressure  (atm) 

A' 

1)  It.Dl.Nj 

water  gas  partial  pressure  in  free 
stream  (atm) 

^lljO.-urf 

water  gas  partial  pressure  adjacent  to 
particle  surface  (atm) 

^ItlKIj.Nj 

^11  ltO,.»iltf 

HB03  partial  pressure  adjacent  to 
particle  surface  (atm) 

boric  oxide  partial  pressure  in  free- 
streatn  (atm) 

heat  release  of  B(s)  -i-  JO;->-  IB.O:,(/) 
(cal/gm-mol) 

heat  release  of  B(0  + J02  -<■  IB2O3(0 
(cal/gm-mol) 

gas  law  constant,  82.06  atm-cm3/ 
gm-mol  0K 

gas  law  constant,  8.316  • 10-  crg/gm- 
mol°K 

molar  rate  of  boron  consumption 
(gm-mol/scc) 

molar  evaporation  rate  of  boric  oxide 
(gm-mol/scc) 

molar  rate  of  removal  of  B.>03  by 
water  reaction  (gm-mol/scc) 
boron  particle  radius  (cm) 
mass  rate  of  generation  of  boric  oxide 
liquid  (gm/scc) 

mass  rate  of  evaporation  of  boric 
oxide  (gm/scc) 

free  stream  gas  temperature  (°K) 
particle  center  temperature  (UK) 
particle  temperature  (°K) 
surroundings  radiation  temperature 
CK) 

particle  surface  temperature  (’K) 
molecular  volume  of  species  j (cm3) 
oxide  layer  thickness  (cm) 
diifusivity  of  gaseous  boric  oxide  in 
nitrogen  (cm-'scc) 

diifusivity  of  1 1 BO_.  in  nitrogen 
(cir.-'secj 

diifusivity  of  water  gas  in  nitrogen 
(cm-/scc) 


35 


256 


MERRILL  K.  KL\0 


A 1 1,,  heat  absorbed  by  reaction  of  ILO 

with  B-Os(/)  (cal-gm-mol) 

A //,„  heal  of  fusion  of  boron  (cal;gm) 

A // v.\ |.  heat  of  vaporization  of  ILO ;,(/) 

(ca!;gm-mol) 

z.  evaporation  coefficient  of  boric  oxide 

liquid 

ct.H  surroundings  absorptivity 

«,  thermal  dilfusivity  of  boron  (cm -/sec) 

<r  Stefan- Boltzmann  constant.  1.354- 

10  1J  cal/cm- sec  K1 
i particle  emissivity 

r Hertz- Knudscn  impingement  factor 

(gm-mol/cm-  atm  sec) 

0 time  (sec) 

pu  boron  density  (gm<cm:‘) 

boric  oxide  density  (gm/cm‘) 

Superscripts 

Overscore  (7*/*,  A ) indicates  quasi-slcady-statc 
value  in  stability  analysis 

Prime  (T'r,  X')  indicates  perturbation  value  in 
stability  analysis 
Subscript 

0 indicates  initial  value 

INTRODUCTION  AND  BACKGROUND 
Recently.  King  (1972)  published  a model  for  the 
ignition  of  single  boron  particles  in  hot  oxygen- 
containing  gas  streams.  This  model  was  used  to 
predict  ignition  times  observed  by  Macck  (1969)  in 
ilnl-Hanie  burner  boron  particle  ignition  studies  and 
to  predict  conditions  required  for  ignition  of  boron 
particles  in  air-augmented  rocket  afterburners. 
Macck  (1969)  observed  in  his  llat-llamc  burner 
studies  that  the  combustion  of  boron  occurs  in  two 
successive  stages.  After  heat  up  to  about  1800- 
2000  K..  the  particle  becomes  luminous,  glows  for 
a short  period  of  time,  fades  out.  and  finally 
rc-ignitcs  to  burn  completely  in  a second  stage 
which  is  much  brighter  and  longer  than  the  first 
stage.  This  author  interprets  the  first  stage  as  an 
ignition  stage  during  which  the  boron  particle  is 
coated  by  a molten  boric  oxide  layer  through  which 
oxygen  must  diffuse  to  react  with  the  boion  to 
provide  reaction  heat  for  vaporization  of  the  oxide 
layer  Before  the  second  stage  of  full-tledged 
combustion  can  occur,  the  oxide  layer  must  be 
completely  removed  from  the  boron  particle. 

The  processes  occurring  during  the  ignition  of  an 
initially  oxidc-co;>ted  boron  particle  in  a hot-gas 
stream  are  complex.  First,  there  is  a heat-up  stage 
during  which  the  comparatively  cold  boron  particle 


is  heated  solely  by  convective  and  or  radiative 
fiux(es)  from  hotter  surroundings.  At  sufficiently 
high  particle  temperature  (approximately  1800  K.) 
non-negligible  self-healing  of  the  particle  by 
exothermic  oxidation  begins  to  be  superimposed 
on  the  convective  and  radiative  fluxes.  Since  the 
particle  is  initially  coated  with  an  oxide  layer,  the 
oxygen  for  this  reaction  must  dilfuse  through  the 
liquid  oxide  layer.  As  the  oxidation  of  the  hoi  in 
core  occurs,  of  course,  it  makes  the  oxide  layer 
thicker,  increasing  diflusionnl  resistance.  At  the 
same  time,  as  long  as  the  particle  temperature 
continues  to  rise,  the  viscosity  of  the  oxide  decreases 
and  consequently  the  dilfusivity  of  oxygen  in  the 
oxide  increases.  In  addition,  boron  oxide  simul- 
taneously evaporates  and  diffuses  away  from  the 
layer  at  a rate  which  depends  on  the  particle 
temperature,  (ending  to  thin  the  layer.  However, 
this  evaporation -is  tin  endothermic  process  which 
tends  to  cool  (or  at  least  lower  the  rate  of  heating) 
the  particle.  Once  the  particle  tcmperatuie  rises 
above  the  surroundings  temperature  tite  convective 
and  radiative  heat  lltixes  turn  negative  and  begin  to 
cool  the  particle.  As  long  as  the  sum  of  the  self- 
heating term  and  convective  radiative  heating 
term  remains  greater  than  the  product  of  the  oxide 
vaporization  rale  and  heat  of  vaporization,  the 
particle  temperature  will  continue  to  rise.  If  this 
situation  persists  to  the  point  where  the  remaining 
oxide  layer  is  sufficiently  thin,  a temperature 
runaway  will  occur,  the  particle  will  limsh  cleansing 
itself  of  oxide,  and  full-fledged  combustion  will 
occur.  If.  on  the  other  hand,  the  sum  of  the  self- 
healing  and  convective  radialive  terms  drops  below 
the  vaporization  heat  demand  beloie  the  surtace  is 
cleaned,  the  particle  will  not  ignite. 

In  King's  (1972)  original  model,  it  was  assumed 
that  oxide  production  is  controlled  hv  diffusion  of 
oxygen  through  the  liquid  oxide  film.  Henry's 
Law.  the  Wilke  correlation  for  diffusivity  of  a gas 
through  a liquid  as  a function  of  viscosity,  and  an 
Arrcnhius  viscosity-temperature  law  weie  used  to 
derive  a two-constant  expression  for  the  oxygen 
delivery  rate  to  the  boron  core,  the  two  constants 
being  evaluated  from  Talley  and  Henderson  (1961) 
data  for  the  oxidation  of  boron  tods.  It  was 
assumed  that  the  removal  of  oxide  from  the  surface 
was  cvaporalion-ralc-controllcd.  with  an  activation 
energy  equal  to  the  oxide  heat  ofevapoi.ition  and  a 
pre-exponential  factor  determined  again  from 
Talley  and  Henderson  data.  Three  nonlinear 
unsteady -state  differential  equations  in  boron  mass, 
oxide  mass,  and  temperature  resulting  from 
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application  of  mass  and  enthalpy  balances  were 
integrated  numerically  for  given  boundary  condi- 
tions and  initial  conditions  to  determine  whether  or 
not  a given  particle  would  ignite  and.  if  so.  how 
long  the  ignition  time  would  be  (both  important 
factors  in  determining  whether  a boron  particle 
will  release  its  heat  in  a finitc-rcsidcncc-timc  after- 
burner). 

The  original  model  of  King  (1972)  has  subse- 
quently been  revised  to  eliminate  several  short- 
comings. Details  of  the  modified  model  development 
and  results  predicted  with  it  are  presented  in  the 
subsequent  sections.  Briefly,  the  changes  arc: 

1 ) In  the  original  model,  the  heat  of  fusion  of 
the  boron,  which-  must  be  supplied  vvhcti  the 
particle  reaches  the  boron  melting  point  of  2450'  K, 
was  not  considered — it  is  included  in  a modified 
treatment  of  the  enthalpy  balance  in  the  new 
model. 

2)  It  was  assumed  in  the  original  model  that 
the  liquic  boric  oxide  removal  rate  was  limited  by 
evaporation  kinetics,  with  Talley  and  Henderson 
data  being  used  toobtain  a pre-exponential  factor 
(and  thus  an  evaporation  coefficient  of  approxi- 
mately 0.006).  This  assumption  has  been  questioned 
by  Mohan  and  Williams  (1972)  in  anothci  boron 
ignition  model  (discussed  later)  in  which  they 
assume  the  oxide  removal  rate  to  be  limited  by 
diffusion  of  the  oxide  gas  away  from  the  particle 
rather  than  by  vaporization  kinetics.  A dimension- 
less group,  A z.  .-".iv  be  examined  to  determine  the 
controlling  nroccss.  If  this  parameter  is  near  unity 
the  resistances  of  the  evaporation  and  diffusion 
steps  arc  nearly  equal  and  neither  resistance  may 
be  neglected,  while  if  it  is  much  larger  than  unity, 
the  process  is  essentially  evaporation  rate  con- 
trolled and  if  it  is  much  less  than  unity,  the  process 
is  essentially  diffusion  controlled.  Substitution  of 
expressions  from  Bird.  Stewart,  and  Lightfoot 
(I960)  for  A and  >•  and  use  of  an  average  gas 
temperature  adjacent  to  the  particle  of  approxi- 
mately 2000  K leads  to: 

k_  _ 1.85  ■ IQ-j  ^ 

at  >•  a Prr 

where  z is  the  boric  oxide  evaporation  coefficient, 
r,.  is  the  particle  radius  in  centimeters,  and  P is  the 
pressure  in  atmospheres.  Parametric  studies 
leading  to  best  fit  of  the  data  of  MaCck  (1969)  as 
discussed  later  indicate  that  a should  be  approxi- 
mately 0.04  this  value  is  consistent  with  a value 
for  boric  oxide  of  0.03  ± 0.01  measured  by 


Bkickburn  (1965).  (The  Blackburn  value  wa* 
brought  to  this  author’s  attention  subsequent  to 
determination  of  the  best  lit  value  of  0.04.) 
Substitution  of  this  value  into  equation  d)  yields 
k'y.v  = 46.2/Pr/,,  where  the  particle  radius  is  now 
expressed  in  microns.  In  Macek’s  experiments,  t lie 
particle  radii  were  17.5  to  22  microns  and  the 
pressure  was  one  atmosphere,  yielding  A yv~~2. 
Under  typical  afterburner  condition*  examined, 
the  particle  radii  ranged  from  I to  10  microns  and 
pressure  ranged  from  1 to  20 atmospheres,  yielding 
A,'xr  values  of  approximately  0.2  to  5.  Accordingly, 
in  the  revised  model  no  limiting  assumption  is 
made  regarding  evaporation  kinetics  or  oxide  gas 
diffusion  control  for  the  oxide  removal  step  -rather. 
thcsc-procc.sscs-.irc  considered-as -scries  resistances. 

3)  Water  appears  to  have  more  beneficial 
effect  on  boron  particle  ignition  than  would  be 
predicted  by  assuming  that  it  diffuses  through  the 
oxide  layer  to  react  with  the  boron  core  in  the 
same  manner  as  oxygen.  Accordingly,  the  moocl 
has  been  modified  to  treat  H.;0  vapor  as  reacting 
by  a difftsion-limned  reaction  with  the  boric  oxide 
[H..0  + RA,(/)->2HBO*l  to  aid  in  removal  of 
the  oxide  layer. 

Mohan  and  Williams  (1972).  have  developed 
similar  (though  not  identical)  equations  for  t lie- 
generation  and  removal  of  a boron  oxide  lu.cr  on  a 
boron  particle  and  the  consequent  variation  in 
particle  temperature  with  time.  However,  t.ither 
than  numerically  integrate  the  resulting  equations 
to  determine  whether  or  not  a particle  ignites  and. 
if  so,  how  long  the  ignition  time  is.  they  have  taken 
a different  approach  to  predict  whether  or  not  the 
particle  will  ignite  for  a given  set  of  boundary 
conditions  (with  no  information  regarding  the 
time  involved).  In  this  technique,  they  convert  the 
boron  and  huric  oxidc  mass  balance  equations  into 
an  equation  for  the  derivative  of  oxide  thickness 
with  time  and  use  the  enthalpy  balance  to  develop 
tin  equation  for  the  derivative  of  particle  tempera- 
ture with  respect  to  time.  They  then  set  the 
derivatives  equal  to  zero  and  solve  the  resulting 
algebraic  equations  for  quasi-stcady-state  values  of 
oxide  thickness  and  particle  temperature  for  ‘he 
given  boundary  conditions  (particle  size,  surround- 
ings temperature,  etc.).  They  next  examine  these 
points  for  stability  using  a linearization  technique 
around  the  quasi-stcady-state  values.  As  part  of 
the  current  v.o-k.  this  author  has  applied  Mohan 
and  Williams’  cchniquc  to  his  equations  for  the 
basic  processes  occurring  during  the  boron  panicle 
ignition.  The  basic  differences  between  the 
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stability  analysts  carried  out  by  mis  author  and  by 
Mohan  and  Williams  arc: 

1)  Monanand  Williams  treat  a planar  geometry 
whereas  this  author  treats  a spherical  geometry. 

2)  Mohan  and  Williams  assume  that  the 
removal  of  bo/ic  oxide  from  the  particle  surface  is 
gas-diffusiomcont rolled  whereas  this  author  treats 
both  evaporation  kinetics  and  gus-dillusion  resist- 
ances in  series. 

3)  Mohan  and  Williams  neglect  heat  loss  from 
the  oxide1  layer  into  the  boron  itself  whereas  this 
author  does  not. 

Results  obtained  using  the  stability  analysis 
technique-  and  using  the  numerical  integration 
treatment  of  the  same  basic  equations  are  compared 
and  discussed  in  a later  section. 

MODF.L  DEVELOPMENT 

The  boron  ignition  model  employed  is  depicted  in 
"figure  I.  If  is  a transient  model  treating  particle 
temperature,  boron  mass,  and  oxide  mass  (and. 
consequently,  oxide  thickness)  as  time-dependent 
parameters.  Boric  oxide  is  generated,  boron  is 
consumed,  and  heat  is  released  by  reaction  of 
oxygen  which  diffuses  across  the  oxide  layer  with 
boron  at  the  boron-boric  oxide  interface.  This 
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HG.  I.  Horon  particle  ignition  model  processes. 


oxide  generation  ,v assumed  to  be  controlled  by  tne 
diffusion  of  the  oxygen  through  the  liquid  oxide 
layer.  (Sample  calculations  indicate  that  in  excess 
of  95  percent  of  the  resistance  of  transport  of 
oxygen  from  the  gas  free-stream  to  the  boron-boric 
oxide  interface  is  in  the  liquid  layer.)  It  is  also 
assumed  that  the  heat  released  In  this  reaction 
(and  in  fact  any  heat  gain.-d  or  lost  by  the  pant-  ic) 
is  distributed  rapidly  throughout  the  particle— that 
is.  it  is  assumed  thai  the  particle  is  isothermal  at 
any  given  time.  An  expression  may  he  easily 
derived  for  the -difference  betwceiuparu.ele  surface 
and  center  temperature  as  a function  of  the  rate  ol 
change  of  surface  temperature: 


Ts  - Trh  = 


(tlTJtlQyr,. 


(2) 


Numerical  calculations  indicate  that  for  horon 
particles  being  ignited  by  hot  "as  streams  except  in 
the  later  stages  of  thermal-runaway,  lor  particles  of 
10  micron  radius.  r/7'„ '<10  < .>  • I()i  K second.  At 
2000  'K.  the  thermal  dilfusivity  of  boron  is 
approximately  5 • 10  •’em-  sec.  Substitution  of 
these  values  into  equation  (2)  yield  .i  -rather 
negligible  temperature  dificrentml  from- surface  to 
center  of  17  K.  Similar  calculations  for  oilier 
particle  sizes  yield  the  same  results  since  the 
maximum  tlT^tlO  value  is  essentially  inversely 
proportional  to  the  square  of  the  particle  radius. 
It  should  be  noted  however  that  this  a,  jmption  of 
tin  isothermal  panicle  is  totally  untenable  for  heat 
fluxes  associated  with  laser-ignition  and  the  model 
therefore  cannot  be  used  for  laser-ignition  studies. 

While  oxygen  is  diffusing  in  through  the  liquid 
oxide  layer  to  convert  boron  to  boric  oxide  and 
: cleasc  heat,  boric  oxide  is  simultaneously  vaporiz- 
ing at  the  outer  surface  and  dilfusmg  out  to  the 
frcc-strcam.  AS  previously  mentioned,  both  the 
finite  kinetics  of  vaporization  and  the  diffusion 
process  arc  considered  as  series  resistances  con- 
trolling this  type  of  oxide  removal.  In  addition, 
in  cases  where  the  hot  gas  stream  contains  water 
gas.  this  iLO  is  assumed  to  remove  oxide  from  the 
outer  surface  by  a diffusion-limited  reaction. 
li..O(g)  + B,Oa(/)  ► 2IIB02(g).  The  model  is 
limited  to  treatment  of  low  1 1,0  concentrations 
(less  than  10  to  20  percent)  through  use  of  the  lean- 
gas  approximation  in  treatment  of  'lie  H.O  and 
IIBOj  diffusion.  Finally,  her l is  added  to  or 
removed  from  the  particle  by  convection  and 
radiation. 

The  careful  observer  may  note  in  examining 
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some  of  the  figures  presented  herein  and  considering 
the  vapor  pressure  versus  temperature  relationship 
of  boric  oxide  that  in  some  cases  the  vapor  pressure 
of  the  oxide  exceeds  the  total  pressure  and  that 
boiling  should  thus  be  considered.  However,  at  the 
times  that  such  situations  occur  the  oxide  layer 
is  only  on  the  order  of  1 00  Angstroms  thick  and-ii 
is  not  clear  to  this  author  that  the  phenomenon  of 
generation  of  unlimited  surface  area  for  vaporiza- 
tion associated  with  boiling  will  occur.  Accord- 
ingly, the  phenomenon  of  boiling  is  not  considered 
and  the  suifacc  area  available  for  transport  of  the 
boric  oxide  from  liquid  to  gas  phase  is  assumed  to 
remain  equal  to  the  particle  geometrical  surface 
area.  For  all  cases  studied,  the  calculated  partial 
pressure  of  boric  oxide  gas  at  'he  surface  does  not 
exceed  approximately  It)  to  20  percent  of  the  total 
pressure  until  well  into  thermal  runaway  since, 
due  to  the  low  evaporation  coefficient  of  the  boric 
oxide,  the  pressure  of  oxide  at  the  surface  remains 
considerably  less  than  the  equilibrium  vapor 
pressure.  (Sec  equation  (3).)  The  dark  period 
observed  by  MaiJck  (1969)  between  the  first  bright 
period  and  the  subsequent  vcrydtrigbt  period  in  the 
history  of  an  igniting  boron  particle  iiriiis  flat-flame 
burner  is  probably  due  to  the  superheated  oxide 
Anally  flashing  off  during  some  part  of  the  thermal 
runaway,  temporarily  blocking  the  inflow  of  oxygen 
to  the  particle  and  causing  a temporary  temperature 
drop- -this  part  of  the  overall  history  is  not  included 
in  the  model. 

With  the  above  as«,,-..|)tiohs  and  limitations, 
the  following-three  equations  (boron  mass  balance, 
oxide  mass  balance,  enthalpy  balance)  may  be 
written  (note  that  three  difleren'  enthalpy  balances 
must  be  employed  depending,  on  whether  the 
particle  temperature  is  less  than,  equal  to,  or 
greater  than  the  melting  point  of  boron.  2450  K): 

d>'/>  _ ftiit.yn  )n  ^ 

(10  4«T)./>„ 

<l±  = (Rnj 2 - Rh  - R„M A/ 
dO  d-rppiij,,, 

dTj.= 

dO 

( R u(Qi:.\)  — I<,.;(A// v,\].)  — R/i  ii  I 
I + 4w(i-/.  4-  A’)'  j 

I x [hCI\  - T,.)  + nrv.,.(T\.  A„  - T).)Y 

( 'i)~r'J.pijCi»its  + 4— I'lT-VpiijOjC piijOj 

(5a) 

(7-V  < 2450  ’«./■=  0) 


if 

JO 


l^lt(C?/,’,v)  — W;.;(A//\  u>)  — R/jSllj/ 

j 4-  4 “(<  /,  4-  A )" 

' X [Ifl  If  //•)  + W?.;;(  r\t  y)i  — 
( ijsrj.pu  A II ), 

( T r = 2450  K.O  < /<  I) 


(5b) 


dTj, 

dO 


( *,AQnxJ  “ R|.;(A //v  ,u> ) - R„  A II „ I 

4-  4n,(/'/.  4-  A-F 

{ x [ht t „ - ?)■)  4-  - T).)V 


)sw-*.p4c,.„i  4*  4nry..\  pii.(,jC|.||s,)( 
(T,.  > 2450  "K./=  I) 


(,5c) 


Equation  5a  expresses  the  rate  of  change  of 
particle  temperature  for  temperatures  below  the 
melting  point  of  boron  (2450  K)  with  df.dO  being 
zero  during  this  period.  Equation  (5b)  expresses 
the  rate  of  change  of  fraction  of  the  boron  melted 
at  the  boron  melting  point.  dT,.l(IO  being  zero  until 
/ reaches  unity.  Finally,  equation  (5c)  expresses 
the  rate  of  change  of  particle  temperature  for 
temperatures  in  excess  of  2450  K.  / being  unity 
and  dfjdO  being  z.ero  in  this  region. 

R,t  is  the  molar  consumption  rate  of  boron. 
R,.:  is  the  molar  rate  of  vaporization  of  boron 
oxide,  and  /?„  is  the  molar  rate  of  removal  of 
boric  oxide  by  the  reaction  ILO  4-  ILO:, (/)— * 
2lil)Os.  (Sec  Nomenclature  for  definition  of 
other  symbols). 

Talley  and  Henderson  (196!)  data  were  used  in 
evaluation  of  the  molar  consumption  rate  of  boron 
(as  limited  by  diffusion  acros:  the  liquid  ILO- 
layer).  With  the  use  of  Henry's  Law.  the  Wilke 
correlation  for  difl'usivity  of  a -gas  through  a liquid 
as  a function  of  viscosity,  and  an  Arrcnhius 
viscosity — temperature  law.  the  following  equation 
form  may  be  derived  for  RK: 


(<) 

Talley  and  Henderson  measured  the  oxide  layer 
thicknesses  at  two  temperatures  in  their  experiments 
conducted  with  boron  rods  at  one  atmosphere 
oxygen  pressure.  In  addition,  they  reported  the 
boron  consumption  rate  per  unit  sample  area  for 
these  two  test  conditions.  These  data  permit 
evaluation  of  the  constants  A',  and  AA  in  equation 
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(Ci).  Tiic  resulting  expression  for  the  molar 
consumption  rate  of  boron  is: 

_ 64,8(1 0~8)(<>  + 


It:  the  development  of  an  expression  for  Rh:  (the 
molar  rate  of  evaporation  of  boric  oxide  from  the 
particle  surface)  both  the  kinetic  limitations  on 
evaporation  rate  and  the  diffusional  limitations  on 
transport  of  boric  oxide  gas  from  adjacent  to  the 
particle  out  to  the  frec-strcam  were  considered  as 
series  resistances — thus, 

// ! ' , jo  ,,  , 

" — y.V\l  ji.Oj  i |t*03,*iirtun7 

4-(r/.  + A)' 

“ kfPijjOj.-iirfjm'  P / ) (^) 

In  this  study,  it  was  assumed  that  the  partial 
pressure  of  boric  oxide  in  the  free-stream  (far  from 
the  particle)  was  zero.  Thus.  thru. manipulation  of 
the  two  equalities  in  equation  (8).  the  expression 
for  may  be  rewritten  as: 

R,:  = 4~(/';.  + XfxrP]h0J(\  + w/*)  (8a) 

Following  substitution  of  expressions  for  the 
Hcrtz-Knudscn  impingement  factor  r,  the  boric 
oxide  vapor  pressure  (/^.Oj).  >'md  k (equations  (9), 
(10).  and  (11))  into  (Sa),  we  arrive  at  equation  (12) 
as  a final  expression  for  RK: 

I / 8 R’T,. 

1 ~ 4RTt.  w(A/U')Bsoj 

= 5.30 T]}  1 mol/cnr  atm  sec  (9) 

P°WH  = 1.5 1 • IOV“"'oft,,'/,‘  6)0) 

_ Q„;o,n-/Nu)  _ l.l8(IQ-')7-1/:-YNu) 

2 RT(r,.  + A')  P(rp  + A') 

1.005  • 10lo(rj>  + 

Rk  ~ 7*/4(1  + (4.50(10 >.P(i>  + A)/7),(Nu))] 

(12) 

Regarding  the  action  of  water  vapor  in  the 
ignition  of  boron  particles,  it  is  assumed  that  the 
water  gas  reacts  in  a dilTusion-limitcd  endothermic 
reaction,  H20(g)  + BaOa(0  — 2HBOs(s>)  at  the 
outer  surface  of  the  oxide  layer,  removing  oxide 
and  absorbing  heat.  Thus,  the  molar  rate  of 
removal  of  oxide  by  this  process,  R„.  may  be 


expressed  as: 

Hii 

4-(/>  + Xf 

~20-„  + X)RTl.{,,hO-"-  ":°-M,r"Uc) 

^nnoj.XifNu) 

— , v.n-r  ('  lllso..»urf.i<<4  I1-5' 

4(r,.  + A )RT,, 

(assuming  that  the  partial  pressure  of  HB02  in  the 
frec-strcam  is  negligibly  small). 

A relationship  between  the  partial  pressures  of 
HaO  and  HBCL  at  the  surface  provided  by  thermo- 
dynamics permits  elimination  of  these  quantities 
from  the  two  cqualitics-of  equation  (13)  to  yield 
with  use  of  equation  (' 4)  (Gilliland  (1934))  for  the 
ditTusi viiies.  equation  (15)  for  R,-' 

n ■'Zt'i 

r\  i n . i A ^ (1 - 


= 4.3-  10 


P(V)*+  V'£) 
**  (A/IK), 


^ (A/IK),  (A/H')st 
R„  = (9.15  • IO“V/,)(Nu)0>  + X)T\? 

x exp  jYs.l  (l  - =p)J  x |^-0. 15 
+J 0.0225  + /W,  exp  [-18.1^1  -~)]] 


The  following  values  were  used  for  other  param- 
eters appearing  in  equations  (3-5),  (7),  (12), 
and  (15): 

M1l'n  = 10.82 

pi,  = 2.33  gniiCtn3 
A/IF„:oa  = 69.64 
Pu.o,  = 1.85  gm 'cm3 
Q it  x ==  146.000  cal/molc 
A/fv.xi*  = 90,000  cal/molc 
XH,,  — 75.000  cal/molc 

It  =--  0.347  • IO  °(Nu)7T;<V(r|.  + .V) 
cal/cin-  sec  K 
s = 0.8 
*/,  = 1.0 

C,.„s  = 0.507  + 7.0  ■ 10  :'Tr  cal/gm  °K 
= 0.438  cal'gm  K. 

XU\,  = 498  cal/gm 
C|>m  = 0.675  cal/gm  K. 

Qitxi—  15 1,000  cal/molc 
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For  all  of  i he  cases  considered  to  date  in  this 
study  the  Nusselt  Number  was  set  equal  to  2. 
That  is.  it  was  assumed  that  the  particle  is  essential- 
ly stationary  with  respect  to  the  surrounding  gases. 
This  is  obviously  not  the  case  for  a particle  released 
into  a high-velocity  gas  stream  for  times  very 
shortly  after  such  release.  For  treatment  of  such  a 
situation,  the  model  may  be  simply  modified  by 
addition  of  a velocity-lag  equation  and  suitable 
Nusselt  Number  versus  Reynolds'  Number  relation- 
ships. For  cases  of  practical  interest,  by  (he  time 
the  particle  is  sufficiently  hot  for  the  reaction  heat 
release  terms  in  equation  (5)  to  be  significant,  the 
particle  velocity  has  equilibrated  sufficiently  with 
the  gas  velocity  to  yield  a Nusselt  Number  of 
approximately  2.  This  analysis  is  meant  to  deal 
only  with  the  ignition  period  (defined  as  the  period 
from  initiation  of  first  particle  glow  to  initiation  of 
the  second  stage  of  combustion),  the  preceding 
heatup  time  having  been  treated  in  the  study  by 
Maceh  ( 1 969).  As  mentioned  earlier,  the  evapora- 
tion coefficient,  x.  was  varied  parametrically  in  the 
early  phases  of  this  study,  a value  finally  being 
chosen  (0.04)  which  gave  good  agreement  between 
prediction  and  Macck's  experimental  data. 

The  above  set  of  equations  [(3  5),  (7).  (12), 
(I5)|  was  converted  to  difference  form  and  pro- 
grammed for  numerical  solution  on  the  computer. 
The  following  quantities  were  supplied  as  input 
parameters:  pressure,  free-stream  oxygen  partial 
pressure,  free-stream  water  partial  pressure,  ambient 
temperature,  radiation  temperature  of  the  sur- 
roundings. initial  oxide  layer  thickness,  initial 
particle  radius,  and  initial  particle  temperature. 
Outputs  included  particle  radius,  oxide  thickness, 
particle  temperature,  oxide  mass,  boron  mass, 
oxide  generation  rate,  oxide  evaporation  rate,  and 
rate  of  removal  of  oxide  by  water  vapor,  till  as 
functions  of  time.  These  calculations  were  per- 
formed for  sufficient  time  to  permit  conclusion 
that  for  a given  set  of  input  conditions  particle 
ignition  would  not  occur,  or  if  it  did  occur,  to 
quantitate  the  ignition  delay  time. 

As  a second  approach  to  calculating  ignition 
limits  (minimum  required  ambient  temperature  for 
ignition)  the  stability  analysis  technique  used  by 
Williams  (1972)  on  a somewhat  different  set  of 
equations  describing  the  various  processes  occurring 
during  boron  ignition  was  applied  to  the  above 
equations.  Numerical  calculations  indicated  that 
the  total  boron  plus  oxide  mass  changed  only  very 
sliuhtly  during  the  ignition  process  for  particle 
sizes  of  interest  and  That  the  stability  analysis  in 


all  cases  of  practical  interest  should  be  applied  at 
particle  temperatures  well  below  2450 VK  (boron 
melting  point).  Accordingly,  equation  (5a)  was 
used  (rather  than  (5b)  or  (5c))  in  this  analysis, 
and  equation  3 was  not  used.  r,.  and  (;•/.  + A') 
being  set  equal  to  their  initial  values.  In  addition, 
use  of  this  analysis  was  limited  to  no-water  cases, 
eliminating  the  need  for  equation  (15).  Finally, 
the  variation  of  y.v  k with  particle  temperature  was 
neglected  as  a second  order  effect.  With  these 
approximations,  equations  (4).  (5a).  (7),  and  (12) 
may  be  reduced  to: 

(IT/,  _ Hc-"Mn  Tr 

,10  ~ X ~ T'/r 

+ CCr , - 7',.)+  £>(TJ,A„-  Tj.)  (16) 

AY  b'T  ■-•S.COO/7*  rv,  — ll.000.Tl* 

= (17) 

JO  X T\'.- 

wherc  .4.  B , C.  D.  /;,  and  /'are  functions  only  of 
the  boundary  conditions  (pressure,  oxygen  partial 
pressure,  ambient  temperature)  and  the  initial 
particle  size,  independent  of  Tt.  and  A': 


A = /,  (Initial  Particle  Size,  P 0.) 
B = f,  (Initial  Particle  Size,  P) 

C =/3  (Initial  Particle  Size,  T„) 
U = ft  (Initial  Particle  Size) 

F.  =h  (/*«,) 

F = f0  (Initial  Particle  Size.  P) 


For  a given  set  of  values  of  boundary  conditions 
and  a gben  particle  size,  equations  (16)  and  (17) 
may  be  set  equal  to  zero  and  solved  for  quasi-steady 
values  of  particle  temperature  (T,.)  and  oxide 
thickness  (A').  The  stability  of  this  solution  can 
then  be  investigated  by  linearization  about  this 
quasi-stcady-statc.  Such  linearization  results  in: 


— = -AfA"  + NT'p 
JO 

(18) 

— = -.S' A"  + UT\. 
,10 

(19) 

where: 

A/  = AT,,e'z--rmr‘>'IX' 

(20) 

..  Ac-s"-6Wirr(,  . 22.600\ 

W"  .V  1 + T,.  ) 

/Jt>— ii.ooo/t,,/  , 44,0u0\ 

+ T',r  127V  T'r  > 

C - 4D7?. 

(21) 
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S = ETl.e-iiXmrr>-:x- 

Fc  SS.MO.T, 


l .EL^(,+21§00) 

V 7"y» 


(22) 


. f£ 

r;,-  hr,.  rj.  ) 

Differentiation  of  equation  (IS)  and  subsequent 
substitution  of  equation  (19)  yields: 


dlT,>  (IT). 

— r + (S  - jV)  —f  + (U- M - SN)Tr  = 0 (24) 

ill)’  (10 

lor  which  the  general  solution  is: 

T).  = A,  + k,emi0 

- (S  - N)  ± \ (T-  A )-  - MUM  - NS) 


(25) 

Inanimation  of  this  solution  indicates  that  for 
stability,  both  roots  (/»,  and  »>,)  must  be  negative, 
which  in  turn  requires  that: 

S > N (26) 

UM  > NS  (27) 

Numerical  substitution  indicates  that  if  the  second 
inequality  is  satisfied,  then  the  first  one  is  also 
satisfied — thus,  for  UM  > NS  the  solution  is 
stable  and  the  particle  docs  not  ignite,  while  for 
l'M  < NS  the  solution  is  unstable  and  the  particle 
ignites.  Algebraic  manipulation  shows  that  the 
inequality  required  for  stability  (27)  reduces  to: 


(~r  - “JrJlQT'/-  -Tt)+  0(7].  - 7-,'u „)] 
t C + 4 DT],  > 0 (28) 


liquations  (16),  (17),  and  (28)  arc  used  in  the 
following  manner  to  examine  ignition  limits  for 
ignition  of  boron  particles  in  a hot  gas  stream. 
1 irst.  the  surroundings  radiation  temperature,  the 
gas  stream  temperature,  pressure,  and  oxygen 
partial  pressure  arc  chosen,  liquation  16  is  then  set 
equal  to  zero  and  solved  for  the  qu.isi-stcady-state 
temperature  values  for  various  particle  sizes.  In 
this  manner,  curves  of  T versus  particle  radius  arc 
generated  for  various  ambient  (gas  stream) 
temperatures  (as  shown  in  Figure  16  for  P=  1 
atmosphere.  /’,,,  = 0.2  atmosphere.  7',tAI,  = 300 
k ).  Lquation  1 7 (set  equal  to  zero)  is  then  used  to 


solve  for  the  corresponding  .V  values.  liquation  (28) 
is  finally  used  to  determine  the  unstable  and  stable 
portions  of  each  7V  curve— the  dotted  line  in 
Figure  16  indicates  the  boundary  between  the 
stable  and  unstable  solutions  (ignition  versus  no 
ignition).  Results  and  interpretation  of  this  analysis 
arc  presented  in  the  next  section. 

RESULTS  AND  DISCUSSION 

In  Figure  2.  plots  of  particle  temperature,  oxide 
layer  thickness,  fraction  boron  melted,  oxide  mass 
generation  rate,  and  oxide  mass  evaporation  rate 
vetstis  time  are  presented  for  a typical  case  m 
which  the  numerical  analysis  predicts  particle 
ignition.  The  particle  treated  is  20  microns  in 
diameter  with  an  initial  oxide  layer  thickness  of 
0.1  micron  and  an  initial  temperature  of  1800  K. 
The  pressure  is  5 atmospheres,  the  oxygen  mole 
fraction  is  0.2.  and  the  ambient  temperature  and 
surroundings  radiation  temperatures  are  both 
2100  K.  As  may  be  seen,  once  the  particle 
temperature  exceeds  approximately  2000  K..  the 
oxide  evaporation  rate  rises  above  the  oxide 
generation  rate,  resulting  in  thinning  of  the  oxide 
layer.  Particle  temperature  continues  to  increase  as 
the  oxide  layer  thins  until  the  boron  melting  point 
is  reached.  At  this  point  the  oxide  thickness  and 
temperature  remain  constant  while  the  boron 
melts  (this  process  requiring  approximately  1 5 per- 
cent of  the  total  ignition  time).  After  the  boron 
has  melted,  the  particle  temperature  resumes  its 
rise  and  the  oxide  layer  thins  enough  that  the 
self-healing  term  (limited  by  oxygen  diffusion 
through  the  oxide  layer)  becomes  so  large  as  to 
cause  a thermal  runaway.  At  this  point,  the  super 
heated  oxide  probably  flashes  off  of  the  particle, 
temporarily  interrupting  the  influx  of  oxygen  and 
resulting  in  the  dark  period  observed  by  Miivick 
(1969)— this  process,  as  mentioned  earlier,  is  not 
included  in  the  model. 

hi  Figure  3,  similar  results  are  presented  for  a 
ease  identical  except  for  reduction  of  the  ambient 
temperature  and  effective  surroundings  radiation 
temperature  bv  100  K.  to  2000  K.  "The  initial 
stages  of  the  process  arc  quite  similar  however, 
heat  losses  to  the  surroundings  when  the  particle 
temperature  rises  above  the  2000  K.  ambient 
temperature  are  sufficient  to  retard  the  evaporation 
rate  sufficiently  that  it  drops  back  to  the  generation 
rate  with  the  result  that  the  oxide  layer  ceases 
thinning.  In  this  ease  (referred  to  as  a degenerate 
ignition  ease)  a stable  quasi-stcady-statc  is  reached. 


r™r*yy&&z  K-!.‘ . P'&za? 


BOROS'  l- ARTICLE  IGNITION  IN  HOT  GAS  STREAMS 


263 


FIG.  2.  Predicted  lime  dependence  of  iniporiant  variables  for  particle  which  ignites. 


the  temperature  levels  out  to  a steady  value,  and 
ignition  docs  not  occur. 

Figure  4 presents  results  for  a third  similar  ease, 
this  time  with  a further  100  °K  reduction  in  ambient 
and  radiation  temperature  to  1900  °K.  In  this  ease, 
the  oxide  evaporation  rate  never  docs  exceed  the 


generation  rate  and  the  oxide  layer  continuously 
thickens,  no  ignition  occurring. 

Macek  (1969)  studied  the  ignition  and  combustion 
of  35  and  44  micron  diameter  boron  particles  in  a 
flat-flame  burner.  The  partial  pressure  of  oxygen 
in  the  burner  flame  was  approximately  0.25  for  most 


c_ 


FIG.  3.  Predicted  lime  dependence  of  important  variables  for  degenerate  ignition  ease. 
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FIG.  4.  Predicted  time  dependence  of  important 
variables  for  particle  which  does  not  ignite. 


of  the  tests,  with  a total  pressure  of  one  atmosphere. 
Water  gas  was  present  in  several  of  the  ‘ests. 
Single  boron  particles  were  injected  into  the  flat- 
llamc  burner  product  stream  and  their  history  was 
observed  with  time-resolved  photography.  The 
particles  were  observed  to  glow  briefly,  die  out  for 
a short  time,  and  finally  glow  brightly  (full  com- 
bustion) until  consumed.  Heat  and  momentum 
transfer  calculations  indicated  that  the  first  glow 
occurred  at  1900  °K,  a temperature  which  the 
model  indicates  would  be  reached  with  negligible 
contribution  from  the  boron-oxygen  reaction. 
Ignition  time  was  measured  as  the  lime  between  the 
appearance  of  this  first  glow  and  the  appearance  of 
full  combustion. 

These  experimental  results  were  used  to  check 
the  ignition  model.  Since  the  initial  oxide  thickness 
on  the  particles  was  unknown,  calculations  were 
performed  with  initial  thicknesses  of  0.02  to  0.10 
microns  (a  range  likely  to  cmcompass  the  actual 
initial  thickness)— variation  over  this  factor  of  5 
resulted  in  only  about  20%  variation  in  the 
predicted  ignition  time.  Predicted  and  experi- 
mental results  for  a scries  of  tests  in  which  no 
water  was  present  and  the  oxygen  mole  fraction  in 
the  burner  product  gases  was  approximately  0.25 
tire  presented  in  Figure  5,  in  the  form  of  ignition 
time  versus  ambient  gas  temperature  for  each  o(  the 
two  particle  sizes  (17.5  and  22  micron  radius)  for 
several  assumed  values  of  the  boric  oxide  evapora- 
tion coefficient  (a).  As  may  be  seen,  use  of  a = 0.04 
gives  quite  good  agreement  between  experimental 
data  and  prediction.  Accordingly,  a value  of  0.04 
was  used  for  a in  all  further  studies. 


Another  test  of  the  model  is  its  capability  to 
predict  the  minimum  ambient  temperature  required 
for  ignition  of  the  particles,  a value  also  measured 
by  Macek  (1969).  For  a dry  gas  stream  at  one 
atmosphere  with  an  oxygen  partial  pressure  of 
approximately  0.25  atmospheres.  Macek  found  the 
minimum  gas  temperature  requited  for  ignition  to 
be  1960-1990  °K  when  the  diluent  used  was 
nitrogen  and  1930  °K.  when  the  diluent  was  argon. 
This  compares  extremely  well  with  the  model 
which  predicts  a minimum  ignition  temperature 
under  these  conditions  of  1950-2000  °K..  (The 
range  of  temperatures  given,  rather  than  a precise 
value,  results  from  the  fact  that  between  1950 
and  20C0°K,  the  model  predicts  ignition,  but  with 
ignition  times  of  greater  than  30  milliseconds,  the 
maximum  particle  residence  time  in  Macek's 
apparatus.) 

As  indicated  in  the  model  development  section, 
an  attempt  was  made  to  incorporate  the  effects  of 
water  gas  on  boron  particle  ignition  into  this 
model.  In  Table  I.  predicted  and  observed  ignition 
times  for  exact  test  conditions  employed  by  Macek 
are  presented  for  all  of  his -tests,  including  those 
with  water  gas  in  the  flat-flame  burner  product 
stream.  As  may  be  seen,  based  on  limited  evidence, 
it  appears  that  the  model  underpredicts  the  effect 
of  oxygen  partial  pressure  on  ignition  time  relative 
to  the  effect  of  amb::nt  temperature.  More 
important,  the  model  grossly  ovcrpreJicts  the 
effect  of  water  vapor  on  the  ignition  time.  More- 
over, with  a water  mole  fraction  of  approximately 
0.15-0.20,  Macek  observes  a minimum  ambient 


2100  2200  2300  2400  2500  2000  2700  2800  2900  300O 

GASTEMPERATURE(*K) 

FIG.  5.  Predicted  and  observed  ignition  times  for 
boron  particles  studied  experimentally  by  Macek 
(1969). 
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TABLE  I.  COMPARISON  OF  MEASURED  AND  PREDICTED  IGNITION  TIMES  FOR 
EXACT  CONDITIONS  USED  BY  MACEK  (1969)  IN  FLAT-FLAME 
BURNER  BORON  PARTICLE  IGNITION  STUDIES. 

PTOTAL=1  A™ 
a = 0.04 

TRADIATION  = 300°K 

ASSUMED  INITIAL  OXIDE  THICKNESS  = 0.1# 

IGNITION  TIME  imiac) 


p02  (atm) 

ph2o  <atm> 

ToJ’K) 

PREDICTED 

OBSERVED 

PARTICLE  RADIUS  =17.2/1 

0.23 

0 

2280 

5.26 

4.4 

0.20 

0 

2430 

3.85 

4.8 

0.23 

0 

2870 

2.24 

3.4 

0.20 

0 

2400 

4.09 

5.0 

0.08 

0 

2510 

2.82 

* 

0.28 

0 

2490 

3.50 

3.6 

0.37 

0 

2450 

3.72 

2.1 

0.19 

0.16 

2240 

1.49 

4.0 

0.21 

0.16 

2330 

1.70 

3.5 

0.19 

0.19 

2430 

1.20 

3.8 

0.20 

0.21 

2640 

1.00 

2.6 

PARTICLE  RADIUS- 22.1/i 

0.23 

0 

2280 

7.99 

5.5 

0.20 

0 

2430 

5.83 

5.7 

0.23 

0 

2870 

3.40 

5.0 

0.37 

0 

2450 

5.61 

3.3 

C.19 

0.16 

2240 

2.30 

7.2 

0.21 

0.16 

2330 

2.76 

5.8 

0.19 

0.19 

2430 

1.88 

6.1 

0.20 

0.21 

2640 

1.55 

5.6 

•IGNITED,  BUT  IGNITION  TIME  NOT  MEASURED 


temperature  required  for  ignition  of  1870-1890 
CK.:  the  model  predicts  ignition  with  this  amount 
of  water  vapor  at  as  low  as  1600 '’K  (the  lowest 
ambient  temperature  tested  in  the  program). 
Thus,  it  appears  that  the  attempt  to  model  the 
effect  of  water  vapor  on  boron  particle  ignition  by 
assuming  a diffusion-limited  reaction,  H..O(g)  + 
li.,03(l)  2H  BOo,  removing  oxide  from  the  surface 


was  unsuccessful.  A possible  explanation  is  that 
such  a reaction  (which  is  quite  endothermic)  is  in 
actuality  kinetics-limited.  However,  this  author  is 
unaware  of  any  data  on  the  kinetics  of  this  reaction. 
Therefore,  the  attempt  to  mode!  the  effect  of  water 
vapor  on  boron  particle  ignition  has  been  dis- 
continued (at  least  for  the  time  being)  and  the 
remainder  of  the  work  conducted  in  this  study  is 
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limited  to  ignition  of  boron  particles  in  dry  gas 
streams. 

The  model  of  boron  ignition  developed  in  this 
study  (minus  considerations  of  tlte  effects  of  water 
vapor)  involves  seven  independent  variables, 
values  of  which  must  be  input  to  the  resulting 
numerical  computer  program  for  prediction  of 
particle  ignition  time.  These  parameters  are: 
initial  oxide  thickness,  initial  particle  temperature, 
initial  particle  size,  ambient  temperature  (surround- 
ings gas  temperature),  effective  surroundings 
radiation  temperature,  pressure,  and  oxygen  mole 
fraction  (or  partial  pressure).  As  part  of  this 
study,  each  of  the  above  independent  parameters 
(except  initial  particle  temperature)  has  been 

:tematically  varied  to  determine  its  effect  on 
whether  the  particle  will  ignite  and,  if  so.  what  the 
ignition  time  will  be.  Results  of  these  parametric 
studies  arc  presented  in  Figures  6-15  and  Table  IF. 

In  Figure  6,  the  effect  of  assumed  initial  oxide 
thickness  over  a five-fold  variation  on  ignition 
time  is  presented.  \s  may  be  seen,  the  smaller 
particle  ignition  times  are  more  sensitive  to  assumed 
initial  oxide  thickness.  With  a one  micron  radius 
particle,  increasing  the  assumed  initial  oxide 
t hickness  from  0.02  microns  to  0. 1 0 microns  approxi- 
mate!) doubles  the  predicted  ignition  time,  while 
for  15  micron  radius  particles,  similar  variation  in 
assumed  initial  thickness  changes  predicted  ignition 
time  only  20  to  25"„. 

The  effects  of  pressure  and  oxygen  mole  fraction 
in  the  gas  free-stream  on  minimum  gas  temperature 
required  for  particle  ignition  are  shown  in  Figures 
7 and  8.  These  calculations  were  all  performed  for 


FIG.  7.  Minimum  gas  temperature  required  for 
ignition  versus  total  pressure  for  various  oxygen 
mole  fractions. 


particles  of  2 micron  radius  with  an  initial  tempera- 
ture of  1800°K  and  an  initial  oxide  thickness  of 
0.02  microns.  Similar  trends  were  noted  with 
other  particle  sizes  and  initial  oxide  thicknesses. 
As  shown  in  Figure  7,  increased  total  pressure 
raises  the  gas  temperature  required  for  ignition  at 
low  oxygen  mole  fractions,  but  lowers  the  required 
gas  temperature  at  high  oxygen  mole  fractions. 
When  tlte  data  are  cross-plotted  as  shown  in 
Figure  8.  it  is  seen  that  the  gas  temperature 
required  for  ignition  is  independent  of  oxygen 
mole  fraction  at  one  and  two  atmospheres  total 
pressure,  but  shows  increasing  dependence  on 
oxygen  mole  fraction  at  higher  total  pressure, 
required  temperature  decreasing  with  increasing 


FIG.  6.  IHfect  of  assumed  initial  oxide  thickness 
on  predicted  ignition  time. 


FIG.  8.  Minimum  gas  temperature  required  for 
ignition  versus  oxygen  mole  function  for  various 
total  pressures. 
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FIG.  9.  Ignition  lime  versus  total  picsstirc  for 
various  oxygen  mole  fractions  at  To.  = 2200  K. 


oxygen  mole  fraction.  These  fairly  complex 
dependencies  result  from  interaction  of  the  effect  of 
oxygen  partial  pressure  on  the  oxide  generation 
rate  and  the  effect  of  total  pressure  on  the  oxide 
removal  rate  (increased  total  pressure  leading  to 
increased  resistance  to  diffusion  of  gas  away 
from  the  particle).  Similar  effects  on  the  ignition 
time  at  a fixed  gas  temperature  as  functions  of 
total  pressure  and  oxygen  mole  fraction  arc  shown 
in  figures  9 and  10.  At  low  oxidizer  mole  fractions, 
ignition  ume  increases  with  increasing  total 
pressure  while  at  high  oxidizer  mole  fractions,  the 
pressure  dependency  is  reversed.  Similarly,  at  one 
atmosphere  total  pressure,  ignition  time  is  predicted 
to  he  independent  of  oxygen  mole  fraction  while 
at  higher  total  pressures,  ignition  time  decreases 
with  increased  oxygen  mole  fraction,  the  sensitivity 


FIG.  10.  Ignition  time  versus  oxygen  mole  fraction 
for  various  total  pressure  at  T»  = 2200  'K, 


TIME  (msec) 


FIG.  II.  Oxide  evaporation  rale  versus  lime  at 
various  total  pressures  at  fixed  oxygen  mole  fraction. 

increasing  continuously  with  increased  total  pres- 
sure. 

The  cause  of  the  increase  in  ignition  time  with 
increasing  total  pressure  at  fixed  oxygen  mole 
fraciion  for  low  oxygen  mole  fraction  (0.1.  0 2)  is 
shown  in  Figures  II  and  12.  where  the  oxide 
removal  rate  and  particle  temperature  are  plotted 
against  time.  As  indicated  on  the  figures,  the  ratio 
of  oxide  evaporation  rale  to  the  evaporation  rate 
that  would  be  observed  in  vacuo  decreases  with 
increased  pressure,  since  the  Jiffusional  resistance 
to  oxide  gas  removal  from  the  vicinity  of  the 
particle  increases.  Since,  for  fixed  oxygen  mole 
fraciion,  the  oxygen  partial  pressure  increases  with 
increased  total  pressure,  early  in  the  process  the 
rate  of  diffusion  of  oxygen  through  the  oxide  layer 
for  reaction  with  the  boron  to  release  heat  mereascs 
with  increased  total  pressure,  leading  to  a faster 


FIG.  12.  Particle  temperature  versus  lime  at  various 
total  pressures  at  fixed  oxygen  mole  fraciion. 
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temperature  rise  as  shown  in  Figure  12.  However, 
the  increased  totai  pressure  also  decreases  the 
oxide  removal  rate  at  any  given  temperature  with 
the  net  result  that  tne  removal  rate  increases  only 
very  little  with  increased  pressure  during  the  early 
stages,  less  than  the  generation  rate  increase. 
As  a result,  the  oxide  layer  thickness  versus  time 
is  higher  in  the  higher  total  pressure  cases,  eventu- 
ally cuu'mg  a decrease  in  heat  release  rate  despite 
the  higher  oxygen  partial  pressure  so  that  the 
temperature  versus  time  curves  cross  in  Figure  12, 
temperature  rising  more  quickly  with  time  in  the 
low  total  pressure  case  than  in  the  high  total 
pressure  case.  This,  combined  with  the  higher 
diffus'onal  resistance  to  evaporation  in  the  high 
pressure  case  causes  the  oxide  remoxal  rate  to  lag 
well  behind  that  in  the  low  pressure  case  in  the 
later  stages.  The  net  result  is  an  increase  in  ignition 
time  with  increased  total  pressure  at  fixed  oxygen 
mole  fraction,  ambient  temperature,  etc.  At 
higher  mole  fractions  of  oxygen  on  the  other  hand 
(0.5  and  1.0)  the  increased  retardation  cf  oxide 
removal  rate  associated  with  increased  total 
pressure  is  overshadowed  by  the  increased  heat 
release  rate  associated  with  the  higher  oxygen 
partial  pressure,  resulting  in  decreased  ignition 
iime  with  increased  total  pressure  at  fixed  oxygen 
mole  fraction,  ambient  temperature,  etc. 

The  effect  of  particle  radius  on  minimum 
ambient  temperature  required  for  ignition  was  also 
studied.  Results,  presented  in  Figure  13.  indicate  a 
decrease  in  minimum  required  gas  temperature 
with  increasing  particle  size  for  panicles  in  the 
1 to  10  micron  radius  size  range.  (The  upper  end 
of  each  bar  indicates  the  lowest  temperature 
examined  which  resulted  in  ignition  while  the  lower 
end  indicates  the  highest  temperature  in\ estimated 
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FIG.  13.  Hfl'ccl  of  panicle  radius  on  minimum 
ambient  temperature  required  for  ignition. 


which  did  nof  result  in  ignition).  The.dcpendency 
noted  with  these  small  particle  sizes  (1  to  10  microns) 
is  associated  with  the  oxide  removal  rate  being 
mainly  exaporatioa  kinetics  limited  (kjv.r  » 1)  for 
the  smaller  sizes. 

In  Figure  14,  predictions  of  ignition  time  versus 
particle  radius  for  particles  which  arc  predicted  to 
ignite  arc  presented  for  several  gas  temperatures 
ranging  from  2050  °K  to  2200  K.  For  the  higher 
temperatures,  the  ignition  time  increases  inono- 
tonically  with  increasing  particle  radius  due  to  the 
increase  in  particle  mass  to  surface  area  ratio 
with  increasing  radius.  At  the  lower  surroundings 
temperatures,  however,  predicted  ignition  time 
decreases  with  increasing  particle  radius  at  the 
lower  sizes  and  then  increases  with  further  particle 
size  radius.  This  region  of  increasing  ignition 
time  with  decreasing  particle  size  is  associated 
with  the  increase  in  minimum  required  gas  tempera- 
ture for  ignition  with  decreasing  particle  size. 

The  effect  of  ambient  temperature  on  predicted 
ignition  time  (radiation  surroundings  temperature 
being  set  equal  to  ambient  gas  temperature)  is 
shown  in  Figure  15.  Ignition  time  may  be  seen  to 
be  quite  sensitive  to  gas  stream  temperature, 
ignition  time  decreasing  with  increased  gas  tempera- 
ture. This  sensitixily  increases  with  decreasing  gas 
stream  temperature  and  increasing  particle  size. 
Finally,  the  effect  of \arying surroundings  radiation 
temperature  at  fixed  gas  stream  ambient  tempera- 
ture is  shown  in  Table  If.  As  would  be  expected, 
ignition  timedeerea:  es  with  increasing  surroundings 
radiation  temperature,  all  other  parameters  being 
held  constant.  In  fact,  for  a gas  stream  temperature 
of  2050  °K  for  the  conditions  studied,  particle 
ignition  will  not  occur  for  a surroundings  radiation 
temperature  of  300  K.  or  2050  K.  but  will  occur 
for  radiation  temperatures  of  2500  Lfv  or  higher, 
ignition  time  dccicasing  markedly  with  lurthei 
increases  in  radiation  temperature  above  2500  K. 
This,  of  course,  implies  that  conditions  for  ignition 
of  an  optically  thick  boron  cloud  or  a boron 
aluminum  or  borommagnesium  cloud  should  be 
much  less  stringent  than  for  ignition  of  a single 
boron  particle  in  a transparent  gas  stream  looking 
at  relatively  cold  radiation  surroundings. 

The  stability  analysis  described  in  the  previous 
section  for  determination  of  whether  or  not  a 
boron  particle  will  ignite  for  given  boundary 
conditions  was  employed  to  generate  the  information 
presented  in  Figures  1(5  and  17  for  a surroundings 
radiation  temperature  of  300  K.  and  a total 
pressure  ol'.onc  atmosphere  (independent  of  oxygen 
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FIG.  14  Idfest  of  particle  radius  on  ignition  time. 


ignition),  l-'or  ambient  temperatures  of  1940  K 
and  higher  there  were  right-hand  and  left-hand 
branches.  The  dotted  line  drawn  through  these 
curves  separates  the  ignition  region  (above  the 
dotted  line)  from  the  no-igmtion  region  (below 
the  dotted  line).  It  is  of  interest  to  compare  these 
curves  to  similar  curves  generated  by  Williams 
(1972).  l-'or  the  larger  particle  sizes,  they  arc  quite 
-imilar  moi  surprising  since  for  the  larger  sizes  the 
dilltisional  resistance  to  11,0a  removal  dominates, 
as  assumed  l>v  Williams).  However,  at  low 
particle  sizes  the  curves  of  Figure  16  and  those 
presented  by  Williams  differ  markedly  in  that 
Williams  pred-ets  no  left-hand  branches  of  ihe 
solid  curves.  This  dilference  results  from  the  fact 
that  he  assumed  diffusion-limned  If.O,  removal 
even  for  small  particle  sizes  whereas  the  model 
presented  herein,  which  treats  the  evaporation 
kinetics  and  didtisi  >n  process  as  scries  resistances 
indicates  that  for  the  smaller  sizes,  the  evaporation 
kinetics  dominate. 

Interpretation  of  the  results  plotted  in  Figure  16 
is  as  follows.  At  any  ambient  temperature  below 
1940  k.  the  particle  will  not  ignite  unless  the 
particle  tempetature  is  initially  greater  than  or 
equal  to  the  value  given  on  the  upper  curve  for 
that  ambient  temperature  for  the  given  particle 
radius.  Though  this  value  of  initial  particle 
temperature  is  a necessary  condition  for  ignition, 
it  is  not  clear  that  it  is  suliiciciti  one  can  imagine 
i case  where  the  initial  particle  temperatinc  exceeds 
the  necessary  value,  but  the  initial  oxide  thickness 


mole  fraction).  For  each  of  several  ambient  gas 
temperatures  equations  (16)  and  (IT)  (with  the 
derivatives  set  equal  to  zero)  were  used  to  generate 
curves  of  ff.  (quasi-steady  state  temperature) 
versus  particle  radius  (solid  curves  in  Figure  16). 
Equation  (28)  was  then  used  to  determine  the 
"stable"  and  "unstable"  poitions  of  these  eui >cs. 
For  ambient  gas  temperatures  below  1940  K. 
upper  and  iowei  curies  were  found  'or  <.  wh  g.o 
temperature  in  all  such  cases,  the  upper  cutve 
repiescnted  unstable  solutions  (ignition)  while  the 
lower  curve  represented  stable  solutions  (no 


I-IG.  15  I fleet  of  gas  icilipcrati  rc  on  boron  p.irtulc 
■eniuon  lime. 
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TABLE  II.  EFFECT  OF  TRAD  =A  TM  ON  IGNITION  TIME 
fp  = 2g,  X0  = 0.02^,  P02  = 1 ATM,  P - 2 ATM,  TP  0 - 1800°  K 


IGNITION  TIME 


Tm(°K) 

tradiation  (°k) 

(imec) 

2200 

300 

C.288 

2200 

0.268 

2500 

0:260 

3000 

0.232 

3500 

0.201 

2100 

300 

0.690 

2100 

0.557 

2500 

0.475 

3000 

0.370 

3500 

0.286 

2075 

300 

1.901 

2075 

0.910 

2500 

0.658 

3000 

0.451 

3500 

0.322 

2050 

300 

NO 

2050 

NO 

2500 

1.442 

3000 

0.600 

3500 

0.372 

particle 

temperature  drops  suffi 

.ntion. 

i>  so  large  mat 
below  the  necessary  \alite  before  the  oxide  layer 
thins  sufficiently  for  thermal  runaway. 

For  ainbicnt  gas  temperatures  of  I ‘W0  K and 
higher,  there  arc  regions  wiicc  there  is  no  particle 
temperature  at  which  <IT,.tl0  and  d\  tIO  can 
simultaneously  equal  zero.  For  instance,  foi  an 
ambient  gas  temperature  of  2000  K.  this  is  the 
case  for  particles  with  radii  between  3.7  and 
42  microns  Thus  a quasi-stcao^  state  cannot 
exist  at  2000  K gas  temperature  for  particles  in 
this  size  region  and  they  will  ignite  for  any  initial 
conditions.  Outside  this  region  (larger  or  smaller 
radius)  a minimum  requirement  on  .aitial  particle 
temperature  exists  for  ignition.  This  rcquiicmcnt 
is  given  by  the  portions  of  the  2000  "K.  gas  tempera- 
ture solid  curves  lying  above  the  dotted  line. 
(Again,  tins  is  a necessary  but  not  necessarily 


be  pert... 
or  not  ig 
heat  uenci 


The  numerical  analvsh  must 
'or  such  cases  to  determine  whether 
will  occur.  However,  since  the 
. terms  in  equations  (5)  and  (16)  are 
much  more  sensitive  to  particle  temperature  than 
to  oxide  thickness  it  seems  likely  t hat  tin  initial 
temperature  more  than  25-50  k above  the 
"necessary"  value  will  result  in  ignition  for  .irtually 
any  initial  oxide  thickness.) 

As  indicated,  for  each  gas  temperature  above 
1940  K,  there  exists  a range  of  particle  sizes  for 
which  ignition  will  occur  independent  of  initial 
conditions,  this  range  being  given  by  the  right 
most  point  of  the  left  branch  of  the  solid  curve  for 
that  gas  temperature  and  the  left  most  point  of  the 
right  branch.  Thus,  the  data  of  Figure  16  may  be 
used  to  generate  a curve  of  minimum  ambient  gas 
temperature  required  to  guarantee  ignition  (for 
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FIG.  16.  Calculated  particle  temperature— panicle 
taditiS  loci  ai tITplilQ  — 0.  tlxjtlO  = 0 for  vai  ions  gas 
temperatures. 

any  initial  conditions)  versus  particle  radius.  Such 
a curve  i->  presented  for  the  case  of  one  atmosphere 
totiil  pressure  and  300  "K.  surroundings  radiation 
lempeiature  in  Figure  17.  As  may  he  seen,  there 
is  a minimum  in  the  curve  at  it  particle  rad:  is  of 
appiosimately  15  microns.  Above  this  size  the 
required  ambient  temperature  to  guarantee  ignition 
under  all  initial  conditions  rises  with  particle  size. 
This  results  from  the  fact  that  for  large  particles 
\\  here  diffusion  controls  the  oxide  removal  rate,  this 
rate,  which  strongly  influences  the  heat  generation 
term,  is  proportional  to  the  particle  radius,  as  is  the 
rate  of  convective  heat  transport  between  particle 
and  surroundings  while  the  radiation  loss  term  is 
proportional  to  the  square  of  the  radius.  Accord- 
ingly. as  may  be  shown  from  manipulation  of 
equation  (16)  (with  proper  substitution  for  A,  B, 
C.  etc.)  with  the  derivative  set  equal  to  zero, 
increasing  particle  radius  requires  a corresponding 
increase  in  gas  temperature  for  the  equality  to  be 
satisfied,  below  a particle  radius  of  approximately 
15  microns  however,  control  of  oxide  removal  rate 
shills  to  vaporization  kinetics  with  the  result  that 
the  heat  generation  term  is  controlled  by  a param- 
eter which  is  proportional  to  the  square  of  the 
particle  radius.  Moreover,  this  term  is  large 
compared  to  the  radiation  term.  Manipulation  of 
liquation  16  shows  that  in  this  regime,  increasing 
particle  radius  requires  a corresponding  decrease 
in  gas  temperature  for  the  derivative  to  be  zero. 

Also  shown  on  Figure  17  arc  results  of  numerical 
calculations  of  minimum  ambient  temperature 
required  for  ignition  for  four  different  particle 
radii.  (The  upper  end  of  (he  bar  indicates  the 


lowest  temperature  tried  for  which  ignition 
occurred  while  the  lower  end  indicates  the  highest 
tested  temperature  at  which  ignition  did  not  occur.) 
As  may  be  seen,  agreement  between  results  of  the 
stability  analysis  procedure  and  the  numerical 
integration  technique  is  excellent. 

The  shape  of  the  curve  in  Figure  17  leads  to  some 
interesting  consequences  for  particles  above  15  /i 
radius.  As  indicated  previously,  the  stability 
analysis  was  performed  using  the  initial  particle 
radius  in  the  various  equations  (neglecting  change 
in  the  radius  during  the  ignition  process).  Accord- 
ingly, the  numerical  analysis  results  for  the  50  p 
radius  particle  shown  in  Figure  17  were  generated 
with  the  particle  radius  held  constant.  If  on  the 
other  hand,  the  radius  is  allowed  to  vary,  the 
particle,  after  an  initial  transient,  quiekiy  avoir's  a 
quasi-steady  state  in  which  the  temperature  and 
oxide  thickness  remain  nearly  constant  and  the 
particle  radius  slowly  decreases.  Eventually  the 
radius  decreases  far  enough  (say  from  50  microns  to 
30  microns)  that  the  minimum  ambient  temperature 
for  ignition  decreases  correspondingly  (in  accord- 
ance with  the  curve  of  Figure  17,  from 
2030  to  1970  "K)  that  the  particle  may  ignite 
even  though  a particle  with  a radius  held  at  50 
microns  would  not  ignite— however,  the  times 
involved  arc  very  long  (thousands  of  milliseconds) 
compared  to  the  times  associated  with  die  normal 
ignition  process. 

A.s  part  of  this  study,  a fundamental  question 
regarding  applicability  of  the  stability  analysis 
approach  was  examined.  Namely,  this  analysis 
involves  solving  for  a quasi-steady-stalc  tempera- 
ture and  oxide  thickness  for  given  boundary 
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conditions  and  examining  tins  point  for  stability: 
thtxqucstion  arises  as  to  whether  one.  starting  from 
an  arbitrary  point  in  the  particle  temperature — 
oxide  thickness  plane  can  expect  to  follow  a path 
in  this  plane  which  will  take  them  close  enough  to 
the  quasi-steady-state  particle  temperature— oxide 
thickness  point  associated  with  the  given  boundary 
conditions  for  the  linearized  stability  analysis  to 
apply.  It  is  quite  unlikely  that  this  can  be  achieved 
for  all  initial  conditions  for  a given  set  of  boundary 
conditions.  However,  results  from  the  numerical 
analysis  showing  essentially  no  dependence  of 
minimum  ambient  temperature  required  for  ignition 
on  initial  conditions  and  the  good  agreement 
obtained  between  predictions  with  the  two  analyses 
indicate  that  for  a wide  range  of  initial  conditions 
the  path  followed  in  the  particle  temperature-oxide 
thickness  plane  passes  close  enough  to  the  quasi- 
steady-state point(s)  associated  with  the  given 
boundary  conditions  to  permit  application  of  the 
stability  analysis.  The  results  of  Figure  18  in 
which  the  paths  followed  by  the  particles  in  the 
particle  temperature-oxide  thickness  plane  arc 
plotted  lor  various  initial  conditions  for  one  cet  of 
boundary  conditions  show  this  clearly.  Over  a 
wide  range  of  initial  conditions,  the  paths  may  be 
seen  to  pass  quite  close  to  the  upper  or  lower 
“Stability”  points,  with  ignition  occurring  in  some 
cases  and  not  in  others.  Thus,  application  of  a 
two-dimensional  stability  analysis  to  boron  ignition 
appears  quite  valid. 

SUMMARY 

A model  of  the  ignition  of  boron  particles  treating 
the  remosal  of  a boron  oxide  layer  from  the  surface 
of  the  boron  particle  through  numerical  integration 
of  difference  equations  describing  the  generation 
and  removal  rate  of  boric  oxide  lias  been  revised 
to  more  accurately  treat  the  various  processes 
boohed.  This  new  model  yields  predicted  ignition 
times  and  minimum  ambient  temperature  require- 
ments for  particles  studied  by  Maeek  (1969)  in  dry 
gas  streams  which  are  in  good  agreement  with  his 
measured  data.  An  attempt  to  treat  the  effect  of 
water  gas  on  boron  particle  ignition  as  being  a 
dilfiision-limited  reaction  of  the  water  with  the 
boric  oxide  coating  to  form  gaseous  II BO*  was 
unsuccessful.  Effects  of  various  parameters  on 
boron  particle  ignition  were  studied  with  this 
model.  Oxide  thickness  was  found  to  moderately 
alfcc!  predicted  ignition  times,  this  effect  decreasing 
with  increasing  particle  size.  The  effects  of  total 
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FIG.  18.  Typical  particle  temperature-oxide  thick- 
ness paths  followed  by  boron  particles  m oxygen- 
containing  gas  stream. 


pressure  and  oxygen  mole  fraction  on  minimum 
ambient  temperature  required  for  ignition  and  on 
ignition  time  wcie  found  to  be  fairly  complex.  At 
low  oxygen  mole  fractions,  increasing  total  pressure 
leads  to  more  difficult  ignition  while  at  high  oxygen 
mole  fractions,  the  reverse  is  true.  At  low  total 
pressure,  ignition  was  found  to  be  independent  of 
oxygen  mole  fraction:  at  high  pressure,  increasing 
oxygen  mole  fractions  lead  to  easier,  faster 
ignition. 

Initial  particle  radius  effects  arc  also  somewhat 
complex.  For  the  case  of  surrounding  radiation 
temperature  equal  to  the  gas  temperature,  mini- 
mum ambient  gas  temperature  required  for  ignition 
increases  with  decreasing  particle  size  for  particle 
radii  below  10  15  microns.  Ignition  time  generally 
decreases  with  decreasing  particle  size  except 
where  affected  by  the  ignition  limit  dependency  on 
particle  size.  Ignition  time  is  predicted  to  increase 
strongly  with  decreasing  gas  temperature.  The 
same  type  of  dependency  on  surroundings  radiation 
temperature  is  predicted  in  tiic  region  of  ambient 
gas  temperatures  which  arc  marginal  for  ignition. 

The  same  equations  used  in  the  numerical 
analysis  have  also  been  used  in  a stability  analysis 
to  define  conditions  (e.g.  ambient  temperature) 
required  for  ignition.  Excellent  agreement  was 
found  between  results  obtained  with  the  numerical 
analysis  and  the  stability  analysis. 


This  study  was  supported  hr  At'OSR  under 
Contract  No.  F44620--71  -C0124  and  was  monitored 
by  Captain  L.  R.  Lawrence. 


52 


BORON  PARTICLE  IGNITION  IN  HOT  GAS  STREAMS 


REFERENCES 

Bird,  R,  B.,  Stewart,  \V.  E.  and  Lighifoot,  E.  N.  (I960). 
Transport  Phenomena,  John  Wiley  & Sons,  NY,  NY. 

Blackburn,  P.  E.  (1965).  Thermodynamics  of  condensed 
and  vapor  phases  in  Ihe  binary  and  ternary  systems 
of  Be-B-O,  Al-B-O,  Si-O,  A'l-Be-O,  and'  Al-B-F, 
Arthur  D.  Little  Final  Report,  ARPA  Project  Code 
No.  9100.  ARPA  Order  No.  315-62. 

Gilliland,  E.  R.  (1934).  hut  Eng.  Client.  22.  1091. 

King.  M.  K.  (1972).  Boron  ignition  and  combustion  in 
air-augmented  rocket  afterburners.  Combustion  Science 
ami  Technology  5,  155-164. 


Mafek.  A.  and  Semple,  J.  (1969).  Combustion  of  boron 
particles  at  atmospheric  pressure.  Combustion  Science 
anil  Technology  1,  181-191. 

Mohan,  G.  .and  Williams,  F.  A.  (1972).  Ignition  and 
combustion  of  boron  in  O./inert  atmospheres.  AHA  J. 
10,  6,  776-783. 

Talley,  C.  P.  and  Henderson,  U.  V.  (1961).  Combustion  of 
elemental  boron.  Proceedings  of  the  4 th  Meeting — 
JANAF-A RPA-NASA  Tlicrtnochcinirai  Panel,  p.  107. 


Received  April  25,  1973 


VII. 
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tional  Symposium  on  Combustion) 
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ABSTRACT 

A detailed  model  of  boron-oxygen-nitrogen  dust  cloud  flames 
Including  consideration  of  the  details  of  boron  particle  ignition  and  the 
effects  of  oxygen  depletion  has  been  developed  and  used  Cor  prediction  of 
flame  speeds  as  functions  of  numerous  parameters.  Reasonably  good  agreement 
between  measured  flame  speeds  and  those  predicted  by  this  model  has  been 
obtained.  In  addition,  a simplified  closed-form  flame  speed  expression  has 
been  developed  and  the  effects  of  the  various  assumptions  used  in  its  develop- 
ment on  predicted  flame  speeds  have  been  examined.  The  detailed  and  simplified 
models  have  both  been  used  to  study  the  effects  of  initial  temperature,  pressur 
initial  oxygen  mole  fraction,  weight  fraction  particles,  initial  particle  si2e, 
initial  thickness  of  the  oxide  coating  on  the  particles,  radiation  feedback 
from  the  post-flame  zone,  and  Nusselt  Number.  Mechanisms  leading  to  the  pre- 
dicted dependencies  are  discussed. 


r7'ir:«Cv  *’  v-~ 

**gg  * 


INTRODUCTION  AND  BACKGROUND 

Boron  is  a particularly  attractive  fuel  for  air-breathing  missile 
engines  due  to  its  high  volumetric  heat  of  combustion.  Among  its  possible 
places  of  application  are  air-augmented  rockets,  slurry  ramjets,  powder  ramjets, 
and  external  burning  missiles.  In  combustor  design  for  these  applications, 
information  regarding  laminar  flame  speed  of  boron-oxygen-diluent  (e.g. , boron- 
air)  dust  clouds  and  the  dependency  of  this  flame  speed  on  various  parameters 
is  important  for  proper  design  of  flame-holding  regions. 

Considerable  effort  has  been  carried  out  over  the  years  on  develop- 
ment of  models  for  prediction  of  laminar  flame  speeds  for  all-gas  systems: 
reviews  of  these  models  have  been  presented  by  Dugger,  Simon,  and  Gerstein^, 
Williams  Von  Karman  and  Stephenson  and  Taylor  among  others.  The 

first  two  references  contain  excellent  discussions  of  what  Williams  calls  the 
"cold  boundary  difficulty".  The  problem  referred  to  is  that  one  cannot  obtain 
a unique  value  of  flame  speed  if  he  allows  finite  rate  chemical  reaction  at  the 
cold  boundary  of  the  system.  This  problem  is  avoided  in  practice  by  choosing 
a temperature  greater  than  the  initial  temperature  of  the  stream  below  which 
chemical  heat  release  is  set  equal  to  xero  and  allowing  the  region  between  the 
initial  temperature  and  the  chosen  temperature  to  be  heated  by  feedback  from 
downstream  regions  (the  ignition  temperature  concept).  A heat  balance  equation 
of  the  following  type  may  then  be  written  at  the  location  of  this  chosen  tempera- 
ture (X  = 0 in  the  coordinate  system  cl.oscn  for  use  in  this  paper): 


(M/A)  C (T.  - T ) = A 
T p i o dx 


+ R. 


(1) 


K = O 
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= mass  flux 

= average  heat  capacity  of  stream  between  temperatures  Tq 
and  T\ 

= ignition  temperature 
= initial  temperature 
= thermal  conductivity  of  stream 
= temperature  gradient  at  location  where  T = 

pc  = o 

Rfl  = radiation  flux  across  X = 0 boundary 

Various  approaches  to  describing  the  temperature  field  for  X > 0 and  thus  per- 
mitting evaluation  of  the  terms  on  the  right-hand  side  of  the  equation  as 
functions  of  the  mass  flux  and  various  independent  parameters  are  described 
for  gaseous  systems  in  References  1-4:  these  permit  determination  of  M/A 

and  thus  the  flame  speed  for  a given  set  of  conditions.  For  given  conditions, 
only  one  value  of  mass  flux  will  permit  integration  of  the  equations  describing 
the  flame  from  the  known  initial  conditions  to  the  known  final  conditions  (2  - 
boundary  eigen-value  problem). 

Several  models  have  also  been  developed  for  prediction  of  flame  speeds 

in  liquid  sprays  (e.g.  Williams  Sundukov  and  Predvoditelev  ^).  For  the 

limiting  case  of  very  fine  droplets  which  atomize  before  significant  chemical 

reaction  occurs,  simple  modification  of  the  gas-flame  models  is  adequate,  while 

for  very  large  droplets  where  negligible  vaporization  occurs  before  ignition, 

the  modeling  of  the  heat  release  distribution  is  also  relatively  simple.  (Both 

(2) 

of  these  limiting  cases  are  described  well  by  Williams.  ) 

In  addition,  several  models  have  been  developed  for  prediction  of  flame 
speeds  in  coal  dust-air  clouds.  Essenhigh  and  Csaba  ^ developed  a model  in 
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which  they  assumed  that  the  ignition  temperature  and  flame  temperature  were 
essentially  equal  and  that  all  heat  transfer  back  to  the  pre-combustion  region 
was  accomplished  by  radiation  from  their  optically  thick  combustion  zone. 

Bhaduri  and  Bandyopadhyay  ^ extended  this  model  by  including  treatment  of 
finite  rate  heat  release  in  the  cloud  (dropping  the  assumption  that  the  ignition 
and  flame  temperatures  were  equal),  but  they  also  neglected  conductive  heat 
feedback  relative  to  radiative  heat  feedback.  In  addition,  they  made  several 
other  assumptions  and  approximations,  including: 

(1)  Temperature  equilibrium  between  particles  and  gas  at  all  positions. 

(2)  No  change  in  gas  density  or  dust  concentration  through  the  cloud. 

(3)  No  ignition  delay  period  (arbitrarily  chosen  ignition  temperature 
independent  of  all  other  variables. 

(4)  Heat  capacities  independent  of  composition  and  temperature. 

For  boron-air  dust  clouds,  preliminary  analysis  indicates  that  conductive  heat 
feedback  cannot  be  neglected  since  the  zone  in  which  the  boron  is  burning  is 
physically  and  optically  quite  thin.  (For  typical  cases,  the  ratio  of  the  area 
subtended  by  boron  particles  in  the  combustion  zone  to  the  total  geometrical 
area  is  on  the  order  of  0.01.)  The  physical  thinness  of  the  combustion  zone, 
typically  less  than  1 millimeter,  indicated  by  this  analysis  also  results  in  a 
high  value  of  temperature  gradient  and  thus  a large  conductive  flux  relative  to 
the  radiative  flux  from  the  combustion  zone.  Downstream  of  the  combustion  zone, 
all  of  the  products  are  gaseous  with  a resultant  low  post-flame  zone  emissivity 
for  reasonable  effective  post-flame  zone  thicknesses  (discussed  further  in  the 
model  development  section).  In  addition,  studies  by  King  indicate  that 

with  boron  there  is  an  appreciable  particle  ignition  delay  time  which  must  be 
included  in  any  accurate  description  of  boron-air  dust  cloud  flames.  Accord- 
ingly, the  two  models  referred  to  above  are  not  applicable  for  prediction  of 
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boron-air  dust  cloud  flame  speeds. 

(10) 

Cassel,  Gupta,  and  Guruswamy  developed  a closed-form  expression 
for  calculation  of  flame  speeds  in  dust  clouds  which  included  radiative  and 
conductive  heat  feedback  through  use  of  an  equation  like  Equation  1 with  sub- 
stitution of  I (T,,  - T . . . ) / Flame  Zone  Thickness!  for  dT/dxl 

*-  flame  ignition  -J  1 x = o 

and  substitution  of  an  expression  for  the  radiation  term  which  summed  the 

radiation  contribution  of  all  of  the  particles  in  the  flame  zone,  assuming 

them  to  all  be  at  the  final  flame  temperature.  The  flame  zone  thickness  was 

2 

estimated  as  a function  of  the  flame  speed  by  use  of  a d ~ burning  law  to 
obtain  the  particle  burning  time  and  estimation  of  the  ratio  of  densities  of 
the  burned  and  unburned  gases  to  permit  conversion  of  this  time  to  a distance. 
This  model  is  much  more  nearly  applicable  to  the  boron-air  dust  cloud  flame 
speed  calculation  problem  than  the  two  previously  referenced  models  — however, 
it  does  have  several  shortcomings  which  are  eliminated  in  the  model  developed 
in  this  study: 


(1)  Particle  ignition  delay  Lime  is  not  treated  analytically;  instead 
an  ignition  temperature  independent  of  other  parameters  is  assumed. 

(2)  Property  (heat  capacity,  density,  thermal  conductivity)  variations 
with  temperature  and  composition  are  not  treated  in  detail. 

(3)  A linear  temperature  profile  with  respect  to  distance  is  assumed 
between  the  ignition  and  burnout  points.  It  is  easily  shown  that 
a continuous  decrease  of  the  temperature  gradient  accompanies  heat 
release  throughout  the  combustion  zone  — thus  the  initial  gradient 
must  be  higher  than  the  average  gradient  over  the  zone. 

(4)  The  effect  of  oxygen  depletion  on  burning  time  is  not  treated. 
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(5)  Change  in  gas  velocity  resulting  from  conversion  of  solid 

reactants  to  gaseous  products  is  not  considered. 

2 

(6)  A d -lav/  is  used  to  describe  the  particle  combustion.  In  the 
pressure-particle  size  regime  of  interest  in  this  study,  Macek^^ 
indicates  that  a d*-law  (kinetics-limited  combustion  rather  than 
diffusion-limited  combustion)  should  be  used  for  boron  particle 
combustion. 

Without  development  of  a detailed  model  treating  the  problems  listed  as  items 
1-5,  it  is  not  clear  how  severely  these  simplifications  affect  flame  speed 
predictions  for  boron-air  dust  clouds.  (The  effect  of  the  last  item,  if  large 
extrapolations  from  particle  sizes  for  which  burn-time  data  are  available  to 
particle  sizes  of  interest  are  used  is  obvious.)  Accordingly,  a complex  computer 
model  employing  a minimum  of  assumptions  and  approximations  has  been  developed  for 
prediction  of  boron-air  dust  cloud  flame  speeds.  In  addition  a much  simpler 
closed-form  expression  using  the  d'-law  for  particle  burning  but  otherwise  using 
much  the  same  approach  as  that  used  by  Cassel,  Gupta,  and  Guruswamy  has  been 

developed  for  comparison. 

MODEL  DESCRIPTION 

Prior  to  development  of  the  detailed  equations  used  in  this  model,  a 
brief  description  of  the  overall  model  will  be  given.  A sketch  of  the  flame 
structure  with  division  into  various  zones  which  are  treated  b>  different  sets 
of  equations  is  given  in  Figure  1.  It  should  be  noted  that  the  model  as  deve- 
loped is  limited  to  fuel-lean  cases,  though  it  could  be  modified  to  also  handle 
fuel-rich  cases.  Preliminary  estimates  indicate  that  Zones  1 through  3 (from 
Tq  + i to  ^34 > where  5 is  a small  number)  are  quite  short  (on  the  order  of  1 to 
5 millimeters)  compared  to  typical  flame  diameters  for  cases  of  interest— 
accordingly  one- dimensional  analysis  may  be  used  for  treatment,  of  these  zones. 
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However,  as  will  be  discussed  later,  the  fact  that  the  post-flame  zone  is 
all  gas  results  in  effective  optical  depths  of  tens  of  centimeters  being 
required  for  attainment  of  emissivities  of  even  0.05  to  0.10.  Thus,  for 
typical  flame  diameters,  multi-dimensional  analysis  of  Zone  4 must  be  carried 
out  for  each  given  geometry  to  obtain  an  estimate  of  the  effective  emissivity 
of  that  zone  as  seen  by  the  preheat  region.  Instead,  it  was  decided  to  simply 
parametrically  examine  the  effects  of  various  assumed  post-flame  zone  effective 
emissivities  in  this  study. 

Zone  1 of  Figure  1 is  the  preheat  region  in  which  the  gas-particle 
mixture  flowing  from  lef t-to-right  entering  at  velocity  Uq  (at  temperature  Tq) 
is  heated  solely  by  conductive  and  radiative  heat  feedback  from  downstream  from 
Tq  to  1600°K.  (1600°K  is  chosen  as  the  beginning  of  the  ignition  zone  since  it 

is  safely  below  the  temperature  at  which  oxide-coated  boron  begins  to  react 

"k  * 

appreciably.)  In  this  analysis,  a value  of  Uq  (or  /A>  the  mass  flux) 

is  guessed.  An  enthalpy  balance  over  Zone  1 then  permits  calculation  of  the 
temperature  gradient  at  the  boundary  of  Zones  1 and  2 (X  = 0) . Since  oxygen 
depletion  effects  are  being  considered  and  since  due  co  nonlinearities  the 
oxygen  mole  fraction  and  temperature  profiles  are  not  similar,  either  the  oxygen 
mole  fraction  or  its  axial  gradient  must  also  be  assumed  at  X = 0 (double  eigen- 
value problem).  Assuming  either  one  permits  calculation  of  the  other  by  a species 
mass  balance  across  Zone  1.  In  this  model,  the  gradient  of  oxygen  mole  fraction 
is  guessed.  Thus,  with  initial  guesses  of  the  flame  speed  and  this  gradient,  all 
conditions  at  X = 0 are  defined,  permitting  the  beginning  of  integration  of  mass. 


species,  and  enthalpy  conservation  equations  thru  Zone  2. 

In  Zone  2,  the  boron  particle  ignition  processes  taka  place.  The 

processes  occurring  during  the  ignition  of  a particle  of  boron  initially  coated 

* Sample  cases  run  using  1400°K  rather  than  1600°K  as  the  temperature  at  the 
end  of  Zone  1 indicated  no  effect  of  this  change  on  calculated  flame  speeds. 
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with  an  oxide  layer  are  complex.  First,  there  is  a heatup  stage  during  which 
the  comparatively  cold  boron  particle  is  heated  solely  by  convective  and/or 
radiative  flux(es)  from  hotter  surroundings.  -Once  the  particle  temperature  is 
sufficiently  high  (approximately  1800°K)  non-negligible  self-heating  of  the 
particle  by  exothermic  oxidation  begins  to  be  superimposed  on  the  convective  and 
radiative  fluxes.  Since  the  boron  particle  is  covered  with  a molten  boric  oxide 
layer  below  about  2100-2500°K,  the  oxygen  must  diffuse  through  the  liquid  oxide 
layer.  As  this  oxidation  occurs,  of  course,  it  makes  the  oxide  layer  thicker, 
increasing  diffusional  resistance.  At  the  same  time,  as  long  as  the  particle 
temperature  keeps  rising,  the  viscosity  of  the  oxide  decreases  and  consequently 
the  diffusivity  of  oxygen  in  the  oxide  increases.  In  addition,  boron  oxide 
simultaneously  evaporates  from  the  layer  at  a rate  which  depends  on  the  particle 
temperature,  tending  to  thin  the  layer.  However,  this  evaporation  is  an  endo- 
thermic process  tending  to  cool  the  particle  (or  at  least  lower  its  rate  of 
temperature  rise).  Once  the  particle  temperature  rises  above  the  surroundings 
temperature  the  convective  and  radiative  heat  fluxes  to  the  particle  will  turn 
negative,  also  tending  to  cool  the  particle.  As  long  as  the  sum  of  the  self- 
heating term  and  the  convective/radiative  term  remains  greater  than  the  product 
of  the  oxide  vaporization  rate  and  the  heat  of  vaporization,  the  particle  tempera- 
ture will  continue  to  rise.  If  this  situation  persists  to  the  point  where  the 
remaining  oxide  layer  is  sufficiently  thin,  a temperature  runaway  will  occur,  the 
particle  will  finish  cleansing  itself  of  oxide,  and  full-fledgea  combustion  will 
occur. 

In  the  integration  of  mass,  species,  and  enthalpy  conservation 
equations  through  Zone  2,  an  analysis  of  single  particle  boron  ignition 

(9) 

as  a function  of  surroundings  conditions  developed  by  King  is  coupled  with 
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these  equations  to  permit  calculation  of  the  removal  of  the  oxide  from  the 
particle  until  thermal  runaway  occurs.  At  this  point,  all  concentrations, 
temperatures,  and  gradients  of  these  parameters  are  stored  for  the  beginning 
of  integration  of  mass,  species,  and  enthalpy  conservation  equations  in  Zone  3. 

In  this  zone,  combustion  of  the  cleansed  boron  particles  by  a kinetics-limited 
(d*-law)  process  with  Macek  data  being  used  for  rate-constant  determination 

is  treated.  The  equations  describing  the  particle  combustion  rate  as  a function 
of  oxygen  concentration  and  particle  size  are  coupled  with  the  mass,  species,  and 
enthalpy  equations  and  the  integration  is  carried  on  through  Zone  3,  which  ends 
when  the  particle  radius  goes  to  zero  (burnout). 

In  cases  where  the  post-flame  zone  emissivity  is  assumed  to  be  zero, 
the  integration  of  the  equations  ends  at  X = X,^.  For  finite  Zone  4 emissivities 
an  energy  balance  is  applied  over  Zone  4 to  correct  the  temperature  at  X = X^  to 
a final  temperature.  The  final  calculated  temperature  and  oxygen  mole  fraction 
are  now  compared  to  values  calculated  from  stoichiometric  and  overall  enthalpy 
balance  equations.  If  there  is  disagreement,  the  initial  guesses  on  Uq  and  the 
oxygen  mole  fraction  axial  gradient  at  X = 0 are  adjusted  ana  the  procedure  re- 
peated. This  cycle  is  continued  until  satisfactory  matching  of  the  final  boundary 
conditions  occurs. 

DETAILED  MODEL  DEVELOPMENT 

Zone  1 - An  energy  balance  may  be  written  over  Zone  1,  in  which  no  heat 
release  occurs,  as  follows  (see  Nomenclature  for  definition  of  symbols): 


l !!i) 

V A / 


(1600-Tq)  cp>Avg-5l  *160u  dx(x  _ o + qRAD,  IN  ~ qRAD,  OUT  ^ 


where  the  total  mass  flux  (M^/A)  of  particles  plus  gas  is  related  to  the  flame 

speed  (U  ) by: 
o 
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As  an  early  part  of  this  study,  a very  simplified  model  (described  later) 
was  derived  and  used  to  calculate  approximate  flame  speeds  and  characteristic 
zone  dimensions  for  typical  cases  of  interest.  7‘hese  values  were  then  used 
in  calculations  to  determine  reasonable  simplifications  of  the  analysis.  One 
such  calculation  indicated  that  the  particle  and  gas  temperatures  are  essen- 
tially equal  at  any  given  location  within  Zone  1 and  thus  that  separate  energy 
balances  for  gas  and  particles  need  not  be  written  in  this  region.  In  addition, 
it  was  found  that  the  burning  zone  is  so  short  as  to  be  optically  thin  (projected 
area  of  all  particles  in  this  region  being  less  than  1 percent  of  the  total  area) 
with  the  result  that  radiation  feedback  from  this  zone  to  Zone  1 is  negligible 
compared  to  conductive  feedback  (less  than  5 percent)  for  cases  of  interest. 
Therefore,  in  consideration  of  radiation  effects  one  is  left  with  radiative 
exchange  between  the  post-flame  (all  gas)  Zone  4 and  the  preheat  Zone  1.  Calcu- 
lation of  Che  effective  emissivity  of  this  post-flame  zone  in  which  the  only 
emitting  species  are  various  boron  oxide  gases  (mostly  I^O^)  is  Qui-te  geometry 
dependent  — accordingly  it  was  decided  to  treat  this  emissivity,  e^,  paramet- 
rically in  this  study.  However,  it  is  instructive  to  estimate  the  post-flame 

zone  emissivity  as  a function  of  its  effective  thickness.  The  approximate 
* 

relationship  between  this  zone  thickness  and  the  emissivity  is  given  by  (through 

(12) 

analogy  to  1^0  and  CO2  gas  emissivity  data  presented  by  McAdams  ): 


AX, 


2 


(4) 


It  should  be  noted  that  this  is  an  empirical  expression  fitting  emissivity 
data  for  path  length-partial  pressure  products  of  5 to  150  cm-atm  and  should 
not  be  extrapolated  outside  this  regime. 
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where  P is  in  atmospheres  and  AX.  in  centimeters.  For  an  initial  particle- 
to-gas  mass  ratio  (k^)  of  0.1  and  a pressure  of  1 atmosphere,  this  equation 
indicates  that  an  effective  post-flame  zone  thickness  of  approximately  400  cm 
is  required  for  an  emissivity  of  0.1,  70  cm  for  an  emissivity  of  0.05,  and 
20  cm  for  an  emissivity  of  0.02.  For  emissivities  in  this  range,  the  afore- 
mentioned simplified  model  indicates  that  radiation  feedback  from  the  post- 
flame zone  to  the  preheat  zone,  while  not  negligible,  does  not  dominate  the 
conductive  feedback.  Based  on  the  above  discussion,  radiative  heat  transfer 
in  and  out  of  Zone  1 was  treated  as  exchange  between  Zones  1 and  4: 

^RAD,  IN  “ ^RAD,  OUT  = °efaB  [Tf  “ F10  (1600)  ] 

where  F^q  is  a correction  factor  for  Zone  1,  ranging  from  approximately  0.3 

for  T = 300°K  to  0.5  for  T = 1000°K,  to  allow  for  the  fact  that  the  entire 
o o 

zone  is  not  at  1600°K.  Equations  2,  3,  and  5 may  be  combined  to  yield: 


dT 

dx 


x 


= o 


P(MW)  _ (1  + k k )u  (1600-T  )C„  . . 

N gl  r oxf  o o P,Avg,l 


RT  A,.™ 
o loOO 


"VbK4  - F10  t1600’*] 



1600 


(6) 


The  average  heat  capacity  over  Zone  1 is  obtained  straightforwardly  by  inte- 
gration of  the  temperature  — dependent  heat  capacity  functions  for  each  ingre- 
dient from  Tq  to  1600°K.  Thus,  given  a guessed  UQ,  dT/dx  |x  _ c for  the  Sas  anc* 
particles  may  be  obtained  from  Equation  6. 

Due  to  the  nonlinear  nature  of  the  equations  for  mass,  species,  and 
energy  conservation  across  Zones  2 and  3,  the  temperature  and  oxygen  mole 
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fraction  profiles  are  not  similar.  Accordingly,  either  the  oxygen  mole  fraction 
or  its  axial  derivative  at  X = 0 must  also  be  guessed  to  start  the  integration 
through  the  downstream  zones  (double  eigen-value  problem).  In  this  model,  the 
value  of  the  gradient  is  guessed  and  the  mole  fraction  itself  is  then  calcu- 
lated from  a mass  balance  on  00  across  Zone  1: 


D1600  1o 

'0Oj  1 + Uo  1600 


O 


(7) 


Thus,  with  selection  of  a temperature  below  which  chemical  reaction 
is  negligible  as  the  boundary  between  Zones  1 and  2 (preheat  and  ignition  zones)  , 
and  guessed  values  of  the  flame  speed  (Uq)  and  the  oxygen  mole  fraction  axial 
gradient  at  x = o,  we  arrive  at  the  initial  conditions  listed  below  needed  to 
begin  integration  of  mass,  species,  and  enthalpy  conservation  equations  for 
gas  and  particles  through  Zone  2. 


T - * 
g,  x = o p,  X = o 


dT 

f 

dx 


o 


(Selected  as  1600°K,  a temperature 
below  which  reactions  are  negligible) 
(Obtained  from  Equation  6) 


dx 


x 


o 


(Guessed) 


Ya  (Obtained  from  Equation  7) 

°2>  x = ° 

r (Known  initial  boron  particle  radius) 

Pi*  = o 

1,  (Known  initial  oxide  coating  thickness 

x = o 

bn  boron  particle) 

Zone  2 - In  Zone  2,  differential  equations  on  gas  temperature,  particle 
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temperature,  and  oxygen  mole  fraction  derived  from  species,  mass,  and  enthalpy 
conservation  equations  are  coupled  with  differential  equations  describing  the 
processes  taking  place  during  the  cleaning  of  the  boron  oxide  coating  from  the 
surface  of  the  boron  particles  (boron  ignition  process)  and  integrated  from  the 
initial  conditions  described  in  the  previous  paragraph  to  the  point  at  which 
the  oxide  thickness  approaches  zero  and  thermal  runaway  occurs  (end  of  Zone  2, 
the  ignition  zone).  The  equations  used  to  describe  the  processes  occurring 
during  boron  particle  Ignition  were  developed  and  described  in  detail  in  Ref- 
erence 9.  Application  of  a differential  enthalpy  balance  to  the  gas  (assuming 
that  all  chemical  heat  release  occurs  in  the  particles  in  this  zone  and  is  only 
transferred  to  the  gas  by  conduction)  leads  to: 


cTT  /dT  \ \ /m  tV 
—j.  + __£  ^ c 

dx  \dx  /Id x/  \ A / 


S T8/-  4nr  “N  h(T  - T ) = 0 
* \dx  / P P g p 


where  the  number  of  particles  per  unit  volume,  N^,  is  given  by: 

3 P (MW)  _ k. 

N 

P 4HRT  r 

g Pi  b 


In  the  derivation  of  Equation  8,  variation  in  thermal  conductivity  with  tempera- 
ture (and  thus  position)  was  allowed  for,  resulting  in  the  second  term.  The 
change  in  molecular  weight  resulting  from  the  consumption  of  oxygen  and  produc- 
tion of  boric  oxide  was  neglected  in  this  ignition  zone  as  being  second  order, 
as  were  changes  in  the  gas  and  particle  mass  fluxes  accompanying  these  reactions. 
(In  general,  only  a very  small  fraction  of  the  boron  is  consumed  during  the 
ignition  phase.) 

Application  of  a differential  oxygen  mole  balance,  including  consumption 
of  oxygen  by  the  boron  particles,  resulted  in: 
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- (MoT/A)  (dY  /dx)  + (MW)-  d[(DP/RT  ) (dY_  /dx)  J/dx-(3R_N_/4) -(MW)  =0  (10) 

S1  u2  &L  S u2  a 1 gi 

where  is  the  molar  consumption  rate  of  boron  per  particle. 

As  with  thermal  conductivity  in  the  enthalpy  equation,  the  variation  in  dif- 
fusivity  with  temperature  (and  thus,  with  axial  position)  was  included.  R,j  is 
the  molar  consumption  rate  of  boron  per  particle. 

A differential  enthalpy  balance  on  the  particles,  including  all  heat 
effects  associated  with  the  reaction  of  oxygen  which  diffuses  through  the  oxide 
layer  to  react  with  the  boron  and  evaporation  of  the  oxide,  and  also  including 
radiative  gains  from  the  post-flame  zone  and  losses  to  the  preheat  zone  may  be 
written  as  (for  < 2450°K.): 


Cpr(t~^+  °5“rN  I!r  2 (Tf4  “ Tn4  ) ~ 
PB'dx  / f B p P \ £ P / 


2 2 
Oo_  N Rr 
B p p 


- Fi0  (l6004) 


4-  N 4 Hr  h T - T I + N 
P P \ 8 P / P 


|)  + Np  rb(qrx)  re  (ahvap)  j 


- 0 


(11) 


Re  is  the  molar  vaporization  rate  of  boric  oxide  per  particle.  At  2450°IC.  the 

boron  particle  begins  to  melt  and  C_„(dT  / dx)  in  the  first  term  must  be  replaced 

. !S  p 

by  _‘il  (df/dx),-  where  f is  the  mass  fraction  of  the  boron  melted.  When  this 
•n 

fraction  reaches  one,  Equation  11  is  again  employed,  with  substitution  of  a new 
heat  capacity  expression  and  a different  value  for  Q , to  allow  for  the  heat  of 
melting.  In  addition  to  the  two  second  and  one  first  order  differential  equations 
above,  two  additional  first  order  differential  equations  describing  mass  bal- 
ances on  the  boron  and  the  condensed  phase  boron  oxide  are  employed  in  Zone  2: 
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The  molar  rates  of  consumption  of  boron  and  evaporation  of  boric  oxide  (per 
particle)  are  given  by  (See  Reference  9): 


64.8(10_8) 
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The  heat  transfer  coefficient  h,  appearing  in  Equations  8 and  9,  the  thermal 
conductivity  A,  and  the  diffusivity  D,  are  calculated  from: 


, = (■.  (Nu) 

r 2 
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A = 0.694  (10“6)  T 0,8 
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dx 
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(16) 
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(20) 


From  the  Zone  2 initial  conditions,,  calculated  as  described  earlier,  the  two 
second- order  differential  equations  (8,  10)  and  the  three  first  order  differ- 
ential equations  (11-13  ) are  numerically  integrated  with  th«  help  of  the 


68 


' ' *'*-**  * > - .*  . * 


^^:^g^sgaa3*i6*w»»iii.  w» 


ancillary  equations  (9,  14-20)  through  Zone  2 unti,l  the  oxide  thickness 
approaches  zero  and  the  particle  temperature  begins  to  run  away.  At  this*  point. 
Zone  2 calculations  are  terminated  and  values  of  r^,  T^,  dT^/dx,  'Y  , aiK*  c^o->"-*x 
are  saved  as  initial  values  for  Zone  3. 

Zone  3 - In  Zone  3,  the  cleaned  boron  particles  undergo  full-fledged 
combustion,  releasing  heat  and  boric  oxide  gas  and  consuming  particulate  boron 
and  oxygen.  Calculations  by  Williams  indicate  that  boron  particles  maintain 

a constant  particle  temperature  of  about  2500°K  during  full  combustion.  Since 
the  ignition  modeling  studies  of  King  indicated  that  igniting  boron  par- 

ticles were  well  into  thermal  runaway  by  the  time  their  temperature  reached 
2500°K,  the  location  at  which  this  particle  temperature  was  achieved  was  chosen 
as  the  termination  point  of  Zone  2.  It  war*  then  assumed  that  this  .particle 
temperature  remained  constant  throughout  Zone  3,  all  heat  release  going  to  the 
gas  phase.  As  a result  differential  equations  required  for  Zone  3 are  limited 
to  equations  in  particle  radius  (r  ),  oxygen  mole  fraction  (Yq  ),  gas  temperature 

(T  ),  and  gas  mass  flux  (M  /A).  (Recall,  that  due  to  the  small  amount  of  con- 
8 8 

•version  of  boron  to  oxide  in  Zone  2,  changes  in  the  relative  amounts  of  gas  .mu 
particle  mass  flux  and  in  effective  stream  molecular  weight  were  ignored  in  that 
zone  — these  effects  are  no  longer  second  order  in  Zone  3 and  must  be  con- 
sidered. A separate  differential  equation  for  particle  mass  flux  is  not  required 
since  its  change  is  simply  the  negative  of  change  in  gas  mass  flux.)  A mass 
balance  equation  on  boron  results  in  the  following  differential  aquation  in 
gas  mass  flux: 


d(M  /A) 


rb  <™>b  np 


69 


'5^*3 > Y-?  • i'/’  ,fvV  «??  ;^''v;.^»3I^ 


where  N can  no  longer  be  calculated  from  Equation  9 due  to  the  non-negli.gibie 
P 

changes  in  molar  gas  flux  and  gas  molecular  weight  affecting  the  gas  density, 
but  must  now  be  calculated  as: 
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P 


3 P (MW)  k (M  /A) 
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(22) 


Application  of  a species  conservation  equation  to  oxygen  results  in  (with  the 


An  enthalpy  balance  written  over  a differential  element  of  Zone  3 (assuming, 
as  mentioned  earlier,  constant  particle  temperature  in  this  zone)  yields: 


/ 


<^XNo 


(24) 


where  R'  is  the  molar  rate  of  consumption  of  boron  per  particle  and  Q'  is  the 
B KA 

heat  release  associated  with  the  reaction  of  1 mole  of  liquid  boron  at  2500°K 
with  3/4  moles  of  09  gas  at  T to  yield  1/2  mole  of  15  gas  at  T^. 

A mass  balance  on  boron  particle  mass  yields: 
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(11) 

Application  of  boron-particle  burning-rate  data  obtained  by  Macek  ' for 
particle  radius-oxygen  partial  pressure  products  below  5 to  10  atmosphere- 
microns  (range  of  interest  for  this  study)  yields  the  following  expression  for 
the  molar  consumption  rate  of  boron  per  particle  in  full-fledged  combination: 


2.0  Rr  2 pn 
p h B 


(MWB) 


(26) 


Equations  17-20  are  used  for  the  thermal  conductivity  (A)  and  diffusiviLy  (!)), 
while  the  gas  molecular  weight  and  heat  capacity  are  estimated  as: 


1 


where  all  of  the  individual  component  specific  heats  are  curve  fit  versus 
temperature. 

Starting  with  the  values  of  r , T , dT  /dx,  Y,.  , and  dYrt  /dx  stored  at 

P g g °2  °2 

the  end  of  the  Zone  2 calculations  and  the  initial  value  of  M /A  , Equations  21, 
23,  24,  and  25  are  simultaneously  integrated  with  the  use  of  ancillary  equations 
17-20,  22,  26-28  across  Zone  3 until  the  particle  radius  r^  goes  to  aero 
(burnout).  In  cases  where  the  post-flame  none  emissivity  is  aero,  calculation 
stops  at  this  point  (since  it  can  be  shown  in  this  case  that  temperature  is 
constant  in  Zone  4)  and  the  final  values  of  gas  temperature  (T  and  oxygen 
mole  fraction  (Yq^  ^ ^ are  comPare^  wit*1-  those  calculated  from  straightforward 
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stoichiometry  and  thermochemical  calculations.  (The  thermochemlcal  calcu- 
lations are  performed  with  no  allowance  for  species  dissociation  and  the  same 
heat  capacity  fits  as  used  in  the  above-described  procedures  for  internal  con- 
sistency.) If  these  do  not  match,  adjustments  are  made  to  the  guesses  on  flame 
speed  (or  equivalently  total  mass  flux)  and  oxygen  mole  fraction  axial  gradient 
at  X = 0 and  the  entire  procedure  is  repeated.  This  loop  is  continued  until 
satisfactory  final  boundary  condition  matches  are  achieved. 

Zone  4 - For  cases  where  a non-zero  er Lssivity  is  assigned  to  Zone  4, 
calculations  must  be  carried  on  through  this  zone  to  account  for  heat  transferred 
out  of  this  zone  back  to  the  preheat  zone.  Since  earlier-mentioned  estimates 
indicated  that  this  zone  must  be  quite  thick  to  have  any  appreciable  emissivity, 
and  thus  temperature  gradients  through  this  zone  are  very  shallow,  conductive 
transport  is  neglected  in  Zone  4 in  this  analysis.  With  this  simplification, 
it  is  readily  shown  that  the  final  temperature  (which  is  the  one  to  be  checked 
against  the  temperature  obtained  from  thermochemical  calculations)  is  approxi- 
mately related  to  the  temperature  of  the  end  of  Zone  3 (T^^)  by 


T = T 

final  g34 


B T£4  - <1600)4  F10 


i T I 


pg34 


(29) 


It  is  also  easily  shown  that  the  oxygen  mole  fraction  does  not  vary  through 

Zone  4:  thus,  Yn  = Yn  . The  same  type  of  final  boundary  condition 

2, final  U2,34 

checks  as  described  above  for  cases  where  e ^ = 0 are  applied  using  T^  ^ in 
place  of  Tg34- 
SIMPLIFIED  MODEL 

As  mentioned  earlier,  a simplified  model  similar  to  thaL  of  Cassel, 
Das  Gupta,  and  Guruswamy^^  except  for  use  of  a d'-law  for  particle  combustion 
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and  a modified  radiation  treatment  consistent  with  the  complex  model  described 

above  was  also  developed  for  comparison  purposes.  The  simple  model  begins  with 

a combination  of  Equations  2,  3,  and  5,  with  (1  + i k ) replaced  by  unity 
4 4 

and  F,10  (1600)  replaced  by  Tq  , leading  to: 


e co.  rr  4 - t * ) + \ — 

f B \ f X0  / T = TI  dx 


p(mwpilc  (t.  - T ) 
S Oo  V i o/ 


T = TI 


where  T^  is  the  minimum  temperature  required  for  ignition  cf  single  boron 
particles  under  static  conditions,  approximately  1900°K.  The  temperature  gradient 
at  the  ignition  point  is  next  approximated  by: 
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The  boron  particle  uiitning  time  is  then  calculated  from  the  d'-law  expression 
used  in  generation  of  Equation  26  (neglecting  oxygen  depletion)  as: 


Sum  ’ k;-  P V0 


Combining  Equations  30  - 32,  one  obtains  the  following  expression  for  flame 
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which  in  the.  aosence  of  radiation  er  = 0 reduces  to: 
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(3*0 


The  following  values  ver°  used 


for  properties  appealing  in 


these  equations: 


X = 


2.5A- 10' 


eal/om  sec  °K 


afi  = 0.t. 

C = 0.2  c.il/fK  *K 
P 


T.  = iVOO'K 

i 


- 30  gm/gm-mol 

- 0.5  cm/atm  sec 


For  each  set  of  conditions  analyzed  with  the  previously  described  detailed 
model,  Equation  33  or  34  was  used  to  calculate  flame  speed  for  comparison  with 
values  calculated  by  that  model. 


RESULTS  AND  DISCUSSION 


Detailed  predicted  profiles  of  gas  temperature,  particle  temperature, 
oxygen  mole  fraction,  particle  radius,  and  oxide  thickness  through  the  ignition  and 
combustion  zones  are  presented  in  Figure  2 for  a typical  boron-air  dust  cloud  case. 
For  this  case,  the  initial  temperature  was  600°K  , the  pressure  was  5 atmospheres, 
the  initial  mole  fraction  oxygen  in  the  gas  was  0.2  (air),  the  weight  fraction  of 
particles  (1  micron  radius,  .02  micron  oxide  thickness)  was  0.06,  the  Nusselt 
Number  was  2 (pure  conduction),  and  the  post-flame  emissivity  was  0.  As  may  be 
seen,  the  oxide  coating  initially  thickens,  due  to  the  rate  of  the  boron  reaction 
with  oxygen  diffusing  through  the  oxide  layer  exceeding  the  evaporation  rate,  until 
a particle  temperature  of  approximately  2000°K  is  reached,  at  which  point  the 
evaporation  rate  rises  above  the  oxide  generation  rate.  Up  to  about  2100°K,  both 
rates  are  sufficiently  slow  that  the  rate  of  heat  transfer  between  particles  and 
gas  is  sufficient  to  keep  both  at  essentially  the  same  temperature.  At  this  point, 
the  oxide  layer  has  been  sufficiently  thinned  and  the  rates  of  oxide  generation 
and  exaporation  raised  to  the  point  where  the  particle  temperature  rises  above  the 
gas  temperature,  eventually  reaching  the  melting  point,  and,  soon  after  melting  of 
the  boron  particle  is  complete,  running  away  (as  the  oxide  layer  thickness  goes 
rapidly  to  zero).  Subsequently,  the  particle  burns  freely  with  the  radius  decreasing 
approximately  linerly  with  distance  until  it  has  been  consumed.  Examination  of  the 
oxygen  mole  fraction  curve  reveals  that  diffusion  from  upstream  of  the  calculation 
starting  point  (Zone  1)  through  Zone  2 results  in  the  oxygen  mole  fraction  through- 
out the  full-fledged  combustion  region  being  nearly  equal  to  the  final  stoichiometric 
value,  indicating  the  importance  of  consideration  of  oxygen  depletion  in  the  modeling 
of  boron-oxygen-nitrogen  dust  cloud  flames. 
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It  should  be  noted  that  the  temperature  of  the  ga..  at  die  point  at  which 
the  boron  particle  begins  full  combustion  is  not  the  static  single  particle  boron 
ignition  temperature  of  approximately  1900°K  (used  in  the  simplified  model),  but 
approximately  2280°K.  Moreover,  examination  of  Table  II  in  which  results  of  various 
cases  are  presented  reveals  that  this  "ignition"  temperature  varies  with  the  different 
independent  variables,  not  being  one  fixed  value  for  all  cases.  While  in  the 
simplified  model  it  is  assumed  that  full  combustion  begins  at  1900°K,  ±t  may  be  seen 
from  Figure  2 that  for  the  case  presented  there,  the  axial  distance  from  the  l'.?00°K 
point  to  the  end  of  the  ignition  zone  is  essentially  equal  to  the  length  of  the 
combustion  zone  itself  (.012  cm  versus  .011cm).  Finally,  the  non-linearity  of  gas 
temperature  versus  distance  in  the  combustion  zone  should  be  noted,  in  comparison 
to  the  assumed  linear  profile  in  the  simplified  model. 

It  is  naturally  of  interest  that  flame  speeds  predicted  with  the  detailed 
model  be  compared  with  experimental  results.  An  extensive  search  for  boron-oxygen- 
nitrogen  dust  cloud  flame  speed  data  unearthed  but  two  data  points,  with  precise 

input  data  available  for  neither  of  them.  These  two  data  points,  supplied  by  Dr. 

<"14) 

Irv  Liebman  ' , were  obtained  in  tests  at  the  Bureau  of  Mines  in  which  the  speed 

of  propagation  of  a laminar  flame  through  a dust  cloud  blown  through  a one-inch 

diameter  tube  was  measured.  The  particle  loadings,  gas  composition,  initial 

temperature,  and  pressure  for  each  test  were  well-defined.  (See  Table  I for  values.) 

However,  the  boron  particle  size  was  not  well  defined  ("under  ten  microns").  For 

the  calculations,  the  diameter  was  assumed  to  be  approximately  5 microns  in  each  case.  In 

addition,  no  information  was  available  on  initial  oxide  thickness.  In  other 

applications  of  the  boron  ignition  model  employed  within  this  model,  use  of  an 

initial  oxide  thickness  of  0.02  microns  was  generally  found  to  give  good  agreement 

with  data:  accordingly,  this  value  was  used  In  the  detailed  model.  (As  will  be 

shown  later,  predicted  flame  speed  is  not  very  sensitive  to  this  parameter.)  From 

the  geometry  of  the  test  apparatus,  the  post-flame  emiss-ivity  was  estimated  to  be  well 
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under  0.01:  at  this  value,  the  effect  of  radiation  feedback  is  at  best  second 

order 

Measured  and  predated  flame  speeds  for  uhe^e  lw.>  eases  are  given  in 
Table  [.  it  should  be  e:rp!;a  that  tuis  is  not  intends  l i > a test  between  die 

simplified  and  detailed  models,  but  rath  ...  as  a che-'k  to  <•«.«*  wnether  these  models 
produce  reasonable  flame-spue J r redictions.  As  may  be  s.  er.,  the  results  art  quite 
encouraging,  the  measured  flame  speeds  h°ing  approximately  12  cm/ sec  in  each  cose 
compa  -ed  to  predictions  of  11. n and  10. A cm/sec  with  the  detailed  model.  Due  to  the 
uncertainties  of  some  of  the  model  input  parameters,  the  data  points  ao  not  provide 
as  rigorous  testing  of  the  model  as  desired;  however,  they  do  indicate  that  the 
model  is  not  unreasonable. 

The  fiarly  close  agreement  between  flame  speed  predictions  by  the 
detailed  and  simplified  models  seen  in  Table  I,  and  latex  in  presentation  of  parameter 
effects  as  predicted  by  the  u.  . models  calls  for  some  discussion.  There  are  three 
major  potential  error-producing  assumptions  associated  wi.cn  the  simplified  model 
(along  with  other  lesser  problems  such  as  neglect  of  temperature  dependence  of  gas 
properties) : 

1.  Neglect  of  the  effect  or  oxygen  depletion.  This  should  result  in 

the  simplified  model  underpredicting  the  combustion  time  and  thus  over- 
predicting the  flame  speed. 

2.  Neglect  of  ignition  delay  time  and  distance.  This  should  also  result 
in  overpredicting  of  the  flame  speed. 

3.  Assumption  of  a linear  ter  oruture  profile  through  the  combustion  zone. 

As  is  well  known,  iieat  release  in  a flame  dominated  by  conductive  heat 
feedback  is  accompanied  by  a decreasing  temperature  gradient  (negative 
second  derivative).  Thus,  the  actual  temperature  profile  in  the  heat 
release  zone  Is  curved,  concave  downward,  and  therefore,  the  temperature 
gradient  at  the  beginning  of  the  zone  (which  determines  the  heat  feedback 
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co  the  preheat  zone)  is  greater  than  the  average  tradient  across 
the  zone.  This  factor  should  result  in  the  simplified  model  under- 
predicting the  flame  speed. 

For  testing  of  the  relative  effects  of  these  errors  and  the  degree  to 
which  the  third  item  cancels  the  first  two  items,  a third  calculation  was  performed 
for  Casel,  in  which  the  third  assumption  was  relaxed,  but  the  first  two  mr~  .itained . 
This  was  done  by  altering  the  detailed  model  to  ignore  all  reactions  in  Zone  2,- 
scart  Zone  3 at  T - 1900°K,  and  to  hold  oxygen  concentration  cual  to  the  initial 
value  throughout  Zone  3.  Results  of  calculations  with  the  detailed  model,  the 
simplified  model,  and  this  third  "in-between"  model  are  shown  in  Figure  3.  As  may 
be  seen,  elimination  of  the  third  assumption  results  in  a considerably  higher  value 
of  (dT/dx)  at  x = 0 c-oir.pared  to  the  "simple"  model  and  thus  to  a much  higher  predicted 
flame  speed.  Further  elimination  of  the  first  and  second  assumptions  (resulting 
in  the  detailed  model),  however,  more  than  compensates,  yielding  a value  of  (dT/dx)  at 
x = 0 even  lower  than  taat  given  by  the  simplified  model  (and  thus  a lower  flame 
speed).  Comparison  of  the  detailed  model  to  the  "in-between"  model  Indicates  a 
drastic  difference  in  predicted  flame  speeds  (4  to  1)  as  does  comparison  of  Che 
"simple"  model  to  the  "in-between"  model  (2-1/2  to  1) : however,  the  errors  are  in 
opposite  directions  and  largely  cancel, resulting  in  the  simplified  model  giving  a 
predicted  flame  speed  only  50%  greater  than  that  predicted  by  the  detailed  model. 

As  may  be  seen  from  Table  II,  this  excess  varies  with  conditions,  ranging  from 
essentially  0%  to  130%  over  the  range  of  conditions  studied. 

As  part  of  this  study,  the  effects  of  eight  independent  parameters  on 
boron-oxygen-nitregen  dust  cloud  flame  speed  ware  examined.  These  parameters 
were  Nusselt  Number,  Initial  Oxide  Thickness,  Post-Flame  Emissivity,  Weight  Fraction 
of  Particles,  Pressure,  Initial  Gas  Stream  Oxygen  Mole  Fraction,  Initial  Boron 
Particle  Radius,  and  Initial  Temperature. 
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Hie  effect  of  Nusselt  Number  on  flame  speed,  as  predicted  by  both  the 
detailed  and  simplified  models  is  shown  in  Figure  4.  Lower  flame  speeds  were 
predicted  with  the  detailed  model.  In  addition,  while  no  dependency  on  Nusselt 
number  was  predicted  with  the-  simplified  model,  the  detailed  model  indicated  a 
decrease  in  flame  speed  with  increasing  Nusselt  Number.  Examination  of  detailed 
profile  plots  similar  to  Figure  2 indicated  that  this  effect  was  caused  by  the 
higher  Nusselt  Number  rerarding  the  takeoff  of  the  Licle  temperature(divergence 
away  from  the  gas  lemperaLure)  lengthening  the  ignition  region  considerably. 

The  predicted  effect  of  initial  oxide  thickness  on  flame  speed  is  shown 
iu  Figure  5.  Again,  considerably  lower  flame  speeds  are  predicted  with  the  detailed 
model  than  with  the  simplified  model.  With  either  model,  flame  speed  is  predicted 
to  decrease  with  increasing  oxide  thickness.  All  of  the  dependency  shown  in  the 
simplified  model  predictions  and  much  of  that  shown  in  the  detailed  model  pre- 
dictions results  from  a decrease  in  flame  temperature  with  increasing  oxide 
thickness  (less  boron  fuel  for  j fixed  weight  fraction  of  parcicles) . (See  Table  ll 
for  flame  temperatures.)  Further  dependence  on  initial  oxide  thickness  shown  by 
the  detailed  model  results  from  the  ignition  zone  length  required  for  remov.il  of 
the  oxide  layer  increasing  with  the  amount  to  be  removed. 

As  shown  in  Figure  6,  and  as  logically  expected , increased  post-flame  emissivity 
results  in  increased  predicted  flame  speeds  with  either  model.  This  result  basically 
derives  from  lower  conductive  heat  feedback  requirements  accompanying  increased 
radiative  feedback.  It  should  be  recalled  from  earlier  discussion,  however,  the 
for  an  initial  par ticla-to-gas  mass  ratio  of  0.1  and  a pressure  of  1 atmosphere, 
an  effective  post-fl"*'.ie  thickness  of  approximately  80  cm  Is  required  even  for  an 
emissivity  of  0.05,  400  cm  for  0.10.  From  simnle  geometric  view-angle  analysis,  the 
ratio  of  effective  post-flame  thickness  to  diameter  for  a circular  jet  flame  is 
approximately  unity — thus,  very  large  diameter  flames  are  required  for  post-fiame 
emissivity  to  have  appreciable  effect  in  the  absence  of  condensed  phases. 
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Increasing  weight  fraction  cf  particles  (on  the  fuel-lean  side  of 
stoichiometric  to  which  this  model  is  currently  limited)  leads  to  predicted 
increases  in  flame  speed  with  either  model.  (Again  lower  values  are  predicted 
with  the  detailed  model.)  This  effect  is  due  to  increasing  flame  temperature  with 
increased  weight  fraction  of  particles.  (See  Figure  7 .) 

Predictions  of  the  effect  of  pressure  on  flame  speed  are  presented 
in  Figure  8.  The  simplified  model  predicts  no  dependency  of  flame  speed  on 
pressure  since  the  effect  of  pressure  on  pre-neat  zone  heat  demand  rate  (pro- 
portional to  pressure  at  fixed  flame  speed)  is  just  offset  by  the  effect  of 
pressure  on  burn  time  and  thus  on  the  temperature  gradient  in  the  combustion 
zone.  With  the  detailed  model,  on  the  other  hand,  flame  speed  is  predicted  to  decreas 
with  increasing  pressure.  The  reason  for  this  is  that  while  increased  pressure 
lowers  die  particle  burning  time,  it  has  little  effect  on  the  ignition  ti.  and 
thus  the  increased  heat  demand  cf  the  preheat  zone  accompanying  pressure  increase 
is  not  totally  offset  by  effects  on  Zones  2 and  3.  Since  the  ratio  of  mass  flux 
to  flame  speed  is  proportional  to  pressure,  it  may  be  seen  that  the  simplified  model 
predicts  a mass  flux  proportional  to  pressure  while  the  detailed  moael  predicts 
mass  flux  increasing  with  increasing  pressure,  but  with  a lesser  dependency. 

As  shown  in  Figure  9,  both  models  predict  increased  flame  speed  with 
increasing  oxygen  mole  fraction,  the  detailed  model  again  resulting  in  lower 
predicted  flame  speeds  than  the  simplified  model.  The  observed  dependency  results 
from  Increased  oxygen  mole  fractions  speeding  up  both  the  ignition  and  combustion 
processes  (particularly  the  latter  since  it  occurs  mostly  at  or  near  the  final 
oxygen  mole  fraction  and  for  a giver,  percentage  increase  in  initial  oxygen  model 
gractior,  M e percentage  increase  in  final  oxygen  mole  fraction  is  greater). 

The  effect  of  increased  initial  boron  particle  radius  is  to  decrease 
predicted  flame  speed  (with  either  model).  This  results  from  increased  particle 
size  causing  increased  combustion  time  and  combustion  zone  length.  This  effect 


is  shown  in  Figure  10. 
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Finally,  increasing  initial  Temperature  strongly  raises  flame  speed, 
as  predicted  by  either  model.  (See  Figure  11.)  This  results  from  two  factors: 

(a)  increased  initial  temperature  reduces  the  preheat  (and  thus  the  heat  feed- 
back) requirement;  and  (b)  increased  initial  temperature  results  in  higher  flame 
temperature,  intensifying  the  heat  feedback.  It  should  bt  noted  that  the  mass  flux 
to  flame  speed  ratio  is  inversely  proportional  to  Initial  temperature,  however, 
since  as  shown  in  Figure  11,  predicted  flame  speed  increase-  more  rapidly 
percentagewise  than  initial  temperature,  then  the  mass  flux  also  increases  with 
initial  temperature. 

SUMMARY 

A model  for  prediction  of  boron-oxygen-r*icrogen  dust  clouds  treating 
the  boron  particle  ignition  and  combustion  pro.~e-.ses  in  detail  was  developed. 

The  problem  was  formulated  as  a double-eigen  value  problem  with  numerical  integra- 
tion of  mass,  species  and  energy  difference  equations  in  combination  with  difference 
equations  describing  the  ignition  and  combustion  processes  through  ignition  ar.d 
combustion  zones.  This  model  was  found  to  yield  predicted  flame  speeds  in  reasonable 
agreement  with  the  limited  bor^n-oxygen-nitrogen  dust  cloud  flame  speed  data  avaiabic. 

In  addition,  a vastly  simplified  closed-form  flame  speed  express!  n was 
developed.  Predicted  flame  speeds  from  this  model  were  found  to  exceed  those 
predicted  from  the  detailed  model  b>  0 to  130%  over  the  range  of  parameters  examined. 
It  was  shown  that  two  of  the  major  error-producing  approximations  associated  with 
the  simplified  model  tend  to  drive  predicted  flame  speeds  down,  while  a third  tends 
Co  drive  them  up,  leading  to  a large  degree  of  cancellation  of  these  errors.  Over 
Che  range  of  conditions  examined,  it  appears  that  the  simple  closed-form  model  can 
be  used  with  a correction  factor  dependent  on  Nusselt  Number  and  pressure  to  obtain 
reasonable  estimates  of  boren-oxygen-nitrogen  dust  cloud  flame  speeds. 
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Effects  of  various  parameters  on  predicted  flame  speed  were  studied. 
Predicted  flame  speeds  were  found  to  increase  with  increasing  initial  temperature, 
decreasing  initial  particle  radius*  initial  oxygen  mole  fraction,  decreasing  pressure, 
increasing  particle  loading  (on  fuel-lean  side  of  scoichiometric) , increasing  post- 
flame e...issivity , decreasing  initial  oxide  coating  thickness,  and  decreasing 
Nusselt  Number. 
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NOMENCLATURE 


oxide  coating  thickness  (cm) 


C , . 

p.Avg.l 


D 

f 


10 


h 

k. 


' oxf 


K 


heat  capacity  (cal/gm  °K) 

average  heat  capacity  between  initial  temperature  and  ignition 
temperature  (cal/gm  °K) 

average  heat  capacity  in  Zone  1 in  detailed  model  (cal/gm  °K) 

2 

diffusivity  of  oxygen  in  oxygen-nitrogen  mixture  (cm“/sec) 
mass  fraction  of  boron  melted 
correction  factor  defined  in  Equation  5 

2 

heat  transfer  coefficient  for  gas-particle  transfer  (cal-/.cm  sec  ®K) 
initial  mass  of  boron  per  unit  mass  of  gas 

initial  mass  of  dxxde  coating  per  unit  mass  of  gas 

constant  in  burning- rate  expression.  Equation  32  (0.5  cm/atnisec) 
molecular  weight 

• 9 

m/a  mass  flux  (gin/cm“sec) 

. 2 

M^/a  Total  mass  flux,  gas  plus  particles  (gm/cm  sec) 

-3 

N number  of  boron  particles  per  unit  volume  (cm  ) 

Nu  Nusselt  Number 

" Pressure  (atm.; 

2 

q .rad, in  radiation  flux  into  preheat  cone,.  Zone  1 ical/cm“sec) 
q rad, out  radiation  flux  out  of  preheat  zone,.  Zone  1 (cal/cm'"sec) 
boron  particle  radius  (cm) 

R gas  law  constant  (82.057  cm  atm/gm-mol  °K) 

2 

net  radiation  flux  into  preheat  zone  (cal/cm  sec) 

Rg  molar  consumption  race  of  boron  per  particle  in  igniton  zone  (gm-mole/sec) 


“R 

MW 


iwiwk-Ww.! 


*&r  \r  ?;».♦  ****;  '? 


' * ;.--> ^ ^ ?’  ^~X3$.p%* gsvsffif 


“i 

re 

Qrx 


<4x 


T 

T, 


^filial 


T. 

j. 


U 


avg 


23 


molar  consumption  rate  of  boron  per  particle  in  combustion  zone 
(gm-mol/sec) 

molar  vaporization  raue  of  boric  oxide  per  particle  (gm-mol/sec) 

heat  release  by  reaction  at  T - T of  .75  moles  of  oxygen  gas  with  1.0 
moles  of  boron  solid  cr  liquid  toPform  0.5  moles  of  I^O,,  liquid  (cal/gm-mol) 

heat  release  bv  reaction  of  .75  moles  of  oxygen  at  T with  1.0  moles  of 
boron  liquid  at  2500°K.  to  form  0.5  moles  of  I^O^  gas8at  T (cai/gm-mol) 

temperature 

flame  temperature  (°K) 

final  temperature  at  end  of  calculations  from  X=0  through  Zone  4 (°K) 
initial  temperature  (°K) 
ignition  temperature  (°R) 

average  velocity  in  combustion  zone  in  simplified  model  (cm/sec) 
flame  speed  [velocity  at  initial  conditions]  (cm/sec) 
distance  from  ocd  of  Zone  1 to  end  of  Zone  2 (cm) 
distance  from  end  of  Zone  1 to  end  of  Zone  3 (cm) 
mole  of  oxygen/moles  gas 


•(,  final 
2 


WB 

ill 


oxygen  mole  fraction  calculated  at  end  of  integration  through  Z jnes 
2 and  3 

boron  particle  absorptivity 


mass  heat  of  fusion  of  boron  (cal/gm) 


vnp 


AX 


flame 


uX. 


molar  heat  of  vaporization  of  boric  oxide  (cal/'  -—mol) 


combustion  zone  thickness  in  simplified  model  (cm) 


effective  radiation  depth  of  post-flame  zone  (cm) 


effective  post-flame  zone  emissivity 
thermal  conductivity  (cal/ciii  sec  °K) 
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.a~- 


3 

density  (gm/cm  ) 

J_2  9 ^ 

Stefan-Boltzmann  Constant  (1.354.10  cal/cm“sec  °K  ) 
particle  burn  time  in  simplified  model  (sec) 

COMMONLY  USED  SUBSCRIPTS 

conditions  at  and  of  preheat  zone  (Zone  1)  in  detailed  model 
initial  conditions 

conditions  at  end  of  Zone  3 (combustion  zone)  in  detailed  model 
boron  property 
boric  oxide  property 

gas  property 
particle  property 
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TABLE  I.  Comparison  of  Predicted  and  Measured  Laminar  Flame 
Speeds  in  Boron-Oxygen-Nitrogen  Dust  Clouds 


TABLE  II.  Effects  of  Various  Parameters  on  Calculated 

Boron-O'xygen-Nitrogon  Dust  Cloud  Flame  Speeds 
Detailed  and  Simplified  Models. 


Gas  Temp, 
at  which 


Parameter  Being  Varied 
And  Value 

Flame 
Temp. 
1°  KV 

Final 

Oxygen 

Mole 

Fraction 

Ignition 
Step  is 
Completed 
(°K) 

Flame  Speed  (U^) 
cm/sec. 

Detailed  Simplified 
Model  Model 

Weight  Fraction  Particles 

.053 

2225 

.086 

2197 

4.73 

10.1 

.060* 

2358 

.074 

2262 

7.67 

11.9 

.065 

2490 

.061 

2314 

9.38 

13.2 

.070 

2620 

.048 

2343 

10.39 

14.3 

.075 

2748 

.035 

2355 

11.30 

15.4 

.080 

2876 

.021 

2387 

12.30 

16.3 

Post-Flame  Eniissivity 

0.00* 

2358 

.074 

2262 

7.67 

11.9 

0.05 

2358 

.074 

2310 

9.52 

14.1 

0.10 

2358 

.074 

2339 

10.99 

16.5 

0.15 

2358 

.074 

2373 

12.57 

18.9 

Initial  Temperature  (°K) 

500 

2274 

.074 

2220 

4.68 

8.7 

600* 

2358 

.074 

2262 

7.67 

11.9 

800 

2532 

.074 

2335 

14.43 

19.6 

1000 

2709 

.074 

2392 

23.19 

30.1 

Pressure  (atm) 

1 

2358 

.074 

2182 

11.45 

11.9 

3 

2358 

.074 

2234 

9.09 

11.9 

5* 

2358 

.074 

2262 

7.67 

11.9 

10 

2358 

.074 

2302 

5.42 

11.9 

Initial  Oxygen  Mole 
Fraction 

0.16 

2353 

.033 

2225 

6.21 

10.5 

0.20* 

2358 

.074 

2262 

7.67 

il.9 

0.24 

2364 

.115 

2285 

8.10 

13.0 

0.30 

2371 

' .177 

2307 

8.90 

14.6 

0.40 

2384 

. 280 

2336 

9.89 

17.1 

Musselt  Number 

2* 

2358 

.074 

2262 

7.67 

11.9 

5 

2358 

.074 

2312 

6.06 

11.9 

Initial  Oxide  Thickness 
(microns) 

0.01 

2406 

.070 

2285 

8.57 

12.3 

0.02* 

2358 

.074 

2262 

7.67 

11.9 

0.05 

2216 

.084 

2193 

4.22 

10.0 

Initial  Particle  Radius 
(microns) 

1.0* 

2358 

.074 

2262 

7.67 

11.9 

2.0 

2406 

.070 

2273 

6.00 

8.7 

3.0 

2421 

.069 

2287 

4.98 

7.2 

5.0 

2434 

.068 

2305 

3.67 

5.7 

* Base  Case:  Weight  Fraction  Particlc*s  = .060 , Initial  Temp.=600°K,  Initial  Oxygen  Mole 

Fraction  0.20,  Initial  Particle  Radius  =-1  micron*  Initial  Oxide  Thickness-0.02  micron 
Prcssuro-5  aim..  Mussel!:  Number -2,  Post-Flame  Kmissivity-0. 
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(REACTIONS 
NEGLIGIBLE 
TO  LEFT  OF 
THIS  PLANE) 

* 

X 


(END  OF 
OXIDE 
COATING 
REMOVAL) 

* 


(PARTICLE 

BURNOUT) 

l 


T = T 


flame 


ZONE  1 

i (PREHEAT  ZONE) 

GAS 

TEMPERATURE 
T = r600°- K 


T = T„ 


= 0 X = 

X23  X = 

X34 

ZONE  3 

ZONE  2 

(PARTICLE 

ZONE  4 

(IGNITION 

COMBUSTION 

(POST-FLAME 

ZONE) 

ZONE) 

ZONE) 

AXIAL  DISTANCE,  X 


FIGURE  1.  DEFINITION  OF  ZONES  CONSIDERED  IN  DEVELOPMENT  OF 

BORON-OXYGEN-NITROGEN  DUST  CLOUD  FLAME  SPEED  MODEL 
(STATIONARY  FLAME  WITH  DUST  CLOUD  ENTERING  FROM  LEFT  AT  VELOCITY  UQ) 


PROFILES  OF  IMPORTANT  VARIABLES  THROUGH  IGNITION  AND 
COMBUSTION  ZONES  AS  PREDICTED  BY  DETAILED  MODEL  FOR  A TYPICAL  CASE 
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riUURE  2.  AXIAL  DISTANCE  FROM  END  OF  PREHEAT  ZONE  (cm) 
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OXYGEN  MOLE  FRACTION 


TEMPERATURE 


MODEL 

I 

II 

III 


UQ  (cm/sec), 

11.9 

30.2 

7.67 


dT/dxix.i.0  ( K/cm) 
53800 

n sopo 

25900 


FLAME  /-END  OF  COMBUSTION 

TEMPERATURE  F ZONE  IN  MODEL  I 


2200  l 


2000 


1800 


1600 


OF  COMBUSTION 
ZONE  IN  MODEL  II 


dT/dx|x  = 0 = 5800°K/cm,  MODEL  ! 


END  OF  COMBUSTION 
ZONE  IN  MODEL  III 

•END  OF  IGNITION 
ZONE  IN  MODEL  III 


dT/dx|x  = 0 = 1 1 6000  K/cm,  MODEL  II 


IGNITION  TEMPERATURE 
USED  IN  MODELS  I AND  II 


25S00°K/cm,  MODEL  III 

i ~.j- 


_L. 


0.005  0.010  0.015  0.020  0.025 

AXIAL  DISTANCE  (cm) 


0.030 


0.035 


FIGURE  3.  COMPARATIVE  TEMPERATURE  PROFILES  FOR  DETAILED 
MODEL  (III),  SIMPLIFIEDsMODEL  (I)  AND  SIMPLIFIED  MODEL 
WITH  ALLOWANCE  FOR  CURVATURE  IN  TEMPERATURE 
PROFILE  ACCOMPANYING  HEAT  RELEASE  (II)  - DEMONSTRATION 
OF  PARTIAL  CANCELLING  OF  ERRORS  ASSOCIATED  WITH 
VARIOUS  ASSUMPTIONS  IN  SIMPLIFIED  MODEL 


FLAME  SPEED  (cm/sec) 


FIGURE  4.  EFFECT  OF  NUSSELT  NUMBER  ON 
PREDICTED  FLAME  SPEED 
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FIGURE  5.  EFFECT  OF  OXIDE  THICKNESS  ON  PREDICTED  FLAME  SPEED 
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FIGURE  7,  EFFECT  OF  PARTICLE  LOADING 
ON  PREDICTED  FLAME  SPEED 


POST-FLAME  EMISSIVITY  = 0 
INITIAL  OXYGEN  MOLE  FRACTION  = 0.20 
INITIAL  TEMPERATURE  = 600°K 
INITIAL  PARTICLE  RADIUS  = 1 micron 

INITIAL  OXSD 
PRESSURE  = 5 
NUSSELT  NUN 

E THICKNESS  = 
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FIGURE  9.  EFFECT  OF  OXYGEN  MOLE  FRACTION  ON 

PREDICTED  FLAME  SPEED 
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INITIAL  OXYGEN  MOLE  FRACTION  IN  GAS  STREAM 
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